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Greek alphabet 


a A alpha t I iota p P tho 
8B B beta « K_ kappa os = sigma 
+ T gamma A A lambda 7 T tau 
6 A delta uw M mu v YT upsilon 
€ E_ epsilon v N mu o @® phi 
¢ Z zeta if Sy x X chi 
n #H eta o O  omicron ¢ WY psi 
009 © theta x I pi w Q omega 


2 SI units (Systéme International d’Unités) 


The SI system is based on six base-units, of which the following four are used in this 
course: 


length: the metre (abbreviated to m); 
mass; the kilogram (abbreviated to kg); 
time: the second (abbreviated to s); 


temperature: the kelvin (abbreviated to K). 


Of these, only the kelvin may be unfamiliar. To express a temperature in kelvins, add 
273.2 to temperature expressed in degrees Celsius or ‘centigrade’ (so the temperature of 
melting ice is 273.2K and that of boiling water is 373.2 KK). The absolute zero of 
temperature is 0K. 


Units for other quantities are given in terms of the base-units, as in the following 
examples: 


velocity: metres per second (ms~!); 
acceleration: metres per second per second (ms~?); 
angular frequency: radians per second (rads~!). 


Notice that ‘per’ is represented by an index —1, and that the usual laws of indices 
apply, so that ‘metres per second per second’ becomes ms~!s~' = ms~*. 


The following combinations of base-units occur commonly enough to have special 
names: 


force: one newton (N) 


W 


one kilogram-metre 

per second per second (kgms~*); 
energy: one joule (J) = one newton-metre (kgm? s~?); 
power: one watt(W) = one joule per second (kgm? s~3). 


Other units are often expressed in terms of these subsidiary units, e.g. 


stiffness of a spring (see Unit 7): newtons per metre (Nm-?); 
dashpot constant (see Unit 8): newtons per (metre per second) (Nm~'s). 


To avoid very large or very small numbers we also use multiple and fractional units, for 
example the kilometre (km) which is 1000 metres. The most important prefixes for 
forming these are given below. 


giga = 1000000000 = 10° (abbreviated to G) 
mega = 1000000 = 10° (abbreviated to M) 
kilo = 1000 = 10° (abbreviated to k) 
milli = 1/1000 = 10-3 (abbreviated to m) 


micro = 1/1000000 = 10-® (abbreviated to 1) 
nano = 1/1000000000 = 10° (abbreviated to n) 


So, for example, the pressure of the atmosphere, which is about 10° Nm~?, is more 
conveniently written 100kN m~?, i.e. 100 kilonewtons per square metre. 


3 Algebra 


3.1 Polynomials 


A polynomial of degree n is an expression of the form 

4,2" + On12"* + +--+ az + a0, 
where n is a positive integer, x is a variable and ap,a),.-.,a, are constants with 
ay #0, A linear polynomial has n = 1, a quadratic polynomial has n = 2 and a 
cubic polynomial has n = 3. 
The roots of a polynomial p(x) are the solutions of the equation p(x) = 0. The number 
of distinct roots is (at most) equal to the degree of the polynomial. If a polynomial 
p(x) has a root a, then it has a factor x — a. If it has a pair of complex conjugate roots 


a and @, then it has a real quadratic factor (x — a)(z —@) = x? + br +c, where 
6=-a-—@and c=aa. 


To factorize a polynomial is to express it as a product of two or more polynomials of 
lower degree. The factors of the difference of two squares x” — a? are (x — a)(x +a); 
those of the perfect square x? + 2ar + a? are (x + a)(x +a), ie. (2 +a)?. 


Completing the square means writing a quadratic polynomial as a sum or difference 
of perfect squares, as follows: 


ar? + br +c = a(x +b/2a)* + k?, 
where k = \/+(c — b?/4a), the sign being chosen to make k real. 


3.2 Summation notation 


Tf a),@2,...,@» is a sequence of numbers, then 


n 
Ya means 4a; + @2 +-+++@n-1 +@n- 


Further, if m is an integer less than n, then 


n 
Ya means Om +@m41 +++ +Gn—1 +@n- 


i=m 


3.3 Factorials and the binomial theorem 


For any positive integer n, we define n factorial, written n!, by 
ni=1x2x3x...x(n-—1)xn. 


The first few are 1! = 1, 2! 3! = 6, 4! = 24. We also define 0! = 1. 


The binomial theorem states that, for any positive integer n, 


(a +6)" =a" + na" b+ mn gn 25? + aes 
nt n—rpr nad n 
+ Fimo BP 4-4 nab"! +6 
“ nt n—rpe 
= age 
Sri(n-r)! 


In particular, 
(a + b)? = a? + 2ab + b*, 
(a +b)° = a® + 3a7b + 3ab? + 6%, 
(a +b)" = a! + 40% + 6076? + dab’ + bi. 
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3.4 Powers and logarithms 


If a > 0, the power function a* has the properties 
a1, 
a*=1/a", 
a® x a¥ =a, 
a*/a¥=a™¥, 
(a*)¥ = a*¥ = (a")*, 


at/¥ = (al/¥)F = (a?)"¥. 


Ifa>0,a#1 and x > 0, the logarithm function log, x has the properties 
log, 1 = 0, 
log, (2~") = —log, x, 
log, (wy) = log, x + log, y, 
log, (2/y) = log, x — loga y. 
log, (x") = y log, 2. 

The power and logarithm functions are inverse functions. Hence 
a=log,y means y=a* 

and aw=a" means y=log,2. 


3.5 Partial fractions 


This technique is useful for breaking up functions of the form P(x)/Q(x), where P and 
Q are both polynomials, into a sum of parts which can be integrated easily. 
A simple case (see Unit 2) is the function 
_ __ar+f 
f(z) = az? +br+c’ 
where b? —4ac > 0. The polynomial ax? + br +c has two real roots: call them x, and 
2, Then we have 
ar+ 3 
f@)= a(x — a )(a — 22)" 
We assume that numbers N;, and N» exist such that 
ar +p Me Np 
a(v—a)(r—a2) 2-2} if r—22 (e# 21,22), 


so that ax + 3 = aNj(x — r2) + aNo(x — 21). 


Equating coefficients of x and the constant terms on both sides of the equation gives 


a pair of simultaneous equations which can be solved for Ny and Nz. So we have 
ea No 
f(e)=— 


ae 
which can be integrated using the table of standard integrals. 


a) | 2-2" 


The general procedure for partial fractions is as follows. 

1, Let the function to be simplified be f(x) = P(x)/Q(x), where P and Q are 
polynomials with real coefficients. 

2. If P(x) has degree equal to or greater than that of Q(), bring it to the form 


P(x) = Pi(x)Q(z) + R(z), 


where P(x) is a polynomial and R(x) is a polynomial of degree lower than that of 
Q(x). Then 
= R(x) 
$2) = Pile) + Gay 

Factorize Q(x) into real linear or quadratic factors. 
4. For each non-repeated linear factor (x — 2;) in Q(x), write a term Nj/(x—2;) in 

the partial fraction expansion. 
5, For each non-repeated quadratic factor x* + jx +c (corresponding to a pair of 
M+ Lic 


Patera in the partial fraction 
i 


complex conjugate roots in Q(.r)), write a term 


expansion. 
6, For each repeated linear factor (2 — x,)"" in Q(x), write terms 
Nia Ni Bee: Nim 
(w= a) | (w— a)? (e-2j)™" 
7. For each repeated quadratic factor (x? + bx + ¢;), write terms 
Miathiae 2 Mim + Lim, 
(x? + bir +c) (a? + bj + ej)™" 


8. Equate the sum of partial fractions to R(x)/Q(x) and multiply both sides of the 
equation by Q(x). Then equate the coefficients of like powers of x on both sides, to 
obtain a set of simultaneous equations for all the unknown constants 
Ni, No,.++)My,Mo,...,L1,L2,..-. 


3.6 Rules of algebra 


In the following summary of the rules or axioms of algebra, the symbols a, b, c stand 
for arbitrary elements of some set S (which can be the set of, for example, real numbers, 
complex numbers, square matrices, n x m matrices, or vectors) and the symbol * stands 
for the binary operation being considered, which may be addition, multiplication, ete. 


Definitions 

Closure: Ifa and b are elements of S, then so is a +b. 

Associativity: (a*b)*c=a+(b*e). 

Identity: There is an ‘identity’ element e such that a+e =e*a =a. 
(For addition it is called 0; for multiplication, 1 or J.) 

Inverse: To each element a there is a unique inverse element @ 


such that a*a@=a*a=e. 
(For addition a is called —a; for multiplication, a~'.) 
Commutativity: a+*b=b«a. 
Distributivity: —a*(boc) = (a+b) o(a*c). 
(over binary (For example, for addition and multiplication of real numbers 
operation 0) a(b+c) = ab+ac.) 


Summary of rules 


Multiplication 
Addition | Real | Complex | Square | Non-square | Vectors: | Vectors: 
(all cases) | numbers | numbers | matrices | matrices dot cross 
product | product 
1. Closure v oe Vi v (iii) i v 
2. Associativity v v v v v x Ne 
3. Identity v v v v x x x 
4. Inverse v (i) (i) (ii) x x x 
5. Commutativity v v v x oe v x 
6. Distributivity x v v v v v v 
(over addition) 


Notes 
(i) Every element except 0 has an inverse. 


(ii) All non-singular square matrices (i.e. those whose rows are linearly independent) 
have inverses, 


(iii) The product of an n x m and an r x s matrix (in that order) exists if and only if 
m =r, in which case the product is an n x s matrix. 


4 Trigonometry 


4.1 Trigonometric functions and their inverses 


Definition for acute | Definition | Range of values 
angles in terms of | in terms of | for function 
Function triangle above sin and cos 
sin® 2 -1<siné<1 
cos 0 ; -l<cos#<1 
aa é sin@ | ai real values | 
a cos 0 
cot 8 . cos@ | Ail real values 
b sin@ 
h zi 
sec@ = — secO>lor<—1 
a cos @ 
h il 
cosec 0 = — cosec@ > 1 or <—1 
b sin® 


Inverse Range of Definition Range of values* 

function | values for x for inverse function 
aresint | -1<a2<1|=6withsind=2 -$<o<t 

arccos 2 -l<xr<1|=@withcoss@=z |0<0<7 

arctanz | Allreale | =@ with tané=xr =t<0<e 

arccotz | All real x =Owithcot@=2 |0<0<7 

arcsec x |x| >1 = 6 with secO = x 0<0<7 and > <O<r 
arceosece | |x| > 1 =O with cosecd = 2 |-7 <0<Oand0<0< 5 


*Sometimes called the principal value range of that inverse function. 


4.2 Two useful triangles 


From these it can be seen that, for example, sin? = — and tan 


v2 


mia 


as 


cory 


4.3 Trigonometric identities 


Pythagorean relations 


sin? @+ cos? @ = 1 
1+ tan? @ = sec? 4 
1 + cot? @ = cosec? 8 


Addition formulas 
sin(a + 3) = sinacos 3+ cosasin 3 
sin(a — 3) = sinacos$ — cosasin# 
cos(a + 3) = cosacos 3 —sinasin 3 
cos(a — 3) = cosacos 3 + sinasin 3 
tana + tan3 


tan(a + 3) = 


tan(a — §) = 


1 


tanatan 3 


tana —tan8 
1+ tanatan 3 
sinacos G = }sin(a + 8) + }sin(a — 8) 
cos asin § = } sin(a + 3) — }sin(a — 8) 
cos acos 8 = } cos(a + 3) + 4 cos(a — 9) 
sin asin 3 = 4 cos(a — 8) — } cos(a + 3) 


1 Be 


CTE 


10 


Double angle formulas 


sin 2a = 2sinacosa@ 
cos 2a = cos? a — sin? a = 1 — 2sin? a = 2cos?a—1 
2tana 


tan 2a = 


1-tan?a 


sin? a = (1 —cos 2a) 


cos* a = $(1 + cos 2a) 


4.4 Polar coordinates 


The point P whose polar coordinates are (r, 4] has Cartesian coordinates (x,y), where 
a=rcosé, y=rsind. 


‘The value of r is always positive (except at the origin, where it is zero). For a given 
point P the value of @ is not unique: we can add or subtract any integer multiple of 27 
and obtain another value for 8 which describes the same point. The value of 6 
satisfying —1 < 0 < 7 is called the principal value of @, and is given by 


ify >0, 


x 
arccos | 
(Zz + *) 


0= 
~arccos | —=—} ify <0. 
x? + y? 


5 Functions and graphs 


5.1 Functions 


Any expression or formula which involves a variable (say x), and whose value is 
uniquely determined by the value of z, is called a function of x. 


If a variable y is a function of z (i.e. is equal to a function of x) then we call x the 
independent variable and y the dependent variable, and we may write y = y(x). 
Here y(z) stands for the function of « (i.e. for the formula involving =). 


A constant is a variable whose value remains the same whatever value we give to the 
independent variable. 


It is also possible to think of a function in a more abstract way, as a mathematical 
object consisting of three parts: 

(i) the domain, which is the set from which the value of x is taken; 

(ii) the codomain, which is a set in which the value of y will lie; 


(iii) the rule which associates to each x in the domain a unique element in the 
codomain, called the image of x. 


li 


If no domain and codomain are specified, then the domain is taken to be the largest 
possible set of real numbers for which the rule makes sense, and the codomain to be the 
set of all real numbers. 


If f and g are two functions, then their sum is a function f +g defined by 
(f +9)(2) = f(z) +9(z) (for all x). 
Moreover, if A and B are any two numbers, then the function Af + Bg is defined by 
(Af + Bg)(x) = Af(z) + Bg(x) (for all x). 
The function f(9(x)) = f 0 g(x) is called the composite function or composition for 
the functions f and g. 
The graph of a function f is the curve in the (x,y)-plane whose equation is y = f(x). 
A constant function is one that assigns the same image to all the elements in its 


domain. Its graph is a straight line parallel to the x-axis, A special case is the zero 
function, which assigns the image zero to all the elements in its domain. 


A linear function is one having the form ax + 6, where a and 6 are constants. The 
constant a is the slope and the constant b is the y-intercept of the graph. 


A continuous function is one such that a small change in x cannot produce a large 

change in f(a) (ie. the graph of f has no breaks or jumps). 

The limit at a of a function f, written lim f(x), is a number Z such that we can make 
ra 


| f(x) — L| as small as we please by choosing any x sufficiently close to a. For a 
continuous function, lim f(x) = f(a). 
ra 


The limit for large x of a function f, written Jim, f(x), is a number L such that we 
can make | f(x) — L| as small as we please by choosing any x sufficiently large. 


5.2 Graphs of some common functions 


See sketches on page 12. 


5.3 Curve sketching 


Suppose we wish to sketch the graph y = f(x), where f(r) is some given function. A 

possible procedure for doing this is given below. 

1, Check if f(x) is any standard function whose graph you already know or a simple 
combination of such functions. 

2. Determine how y behaves when z is very large and positive and when x is very 
large and negative. 

3. Find where the curve crosses the x- and y-axes, if at all. 

4. Look for any values of x at which f(x) is undefined, and examine the behaviour of 
f(x) near these values of x. 

5. Find if there are any local maxima, local minima or points of inflection in the given 
domain. Find these values of f(x) and the associated values of xr. 

6, Transfer the information found in Steps 3 and 5 to a sketch graph, and use what 
you know about the general shape of the graph to help you join the points with a 
smooth curve, taking into account any information found in Steps 2 and 4. If you 
are still unsure about any parts of the curve, choose suitable values of x and plot 
the corresponding points (., f(.)) before completing the curve. 
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Sketch graphs of some common functions 


aS 48 AS 
ey 
2 Soe 8 
S H 
_ 4 4 
g 
< < 
= = 
44 ib 
au 
Cn (0) cn ae 
—— < 
> = 
Se ee Lo 
Ay 4s 
A g 
mth 
> = 3 
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6 Differentiation 


(Only ordinary differentiation is considered here. For partial differentiation, see the 
summary of Unit 25.) 


6.1 Notation and terminology 


If f is a function, its derived function or derivative f’ is defined by 
Hn) = Yim Lett =f) 
1S a 
The process of calculating f‘(x) from f(:r) is called differentiation of f(x) with 
respect to 2. Differentiation with respect to x is denoted by the symbol ee written to 
the left of the expression or variable which is differentiated, so that, for example, 


H(z) means f(x), 
dx 


and if y is a variable that stands for f(x), then # is a variable that stands for f’(x). 
dy 
ey 
When the independent variable is the letter xr, we often use a prime in place of 4, so 
that 4/ me ay When the independent variable is ¢ (time), we often use a dot in 
place of an” that @ means = 


The derivative of a derivative is called a second derivative. For example, the second 
derivative of the function f, denoted by f”, is the (first) derivative of f’, defined by 


f(e+h)= fz) 


To save space, we often print dy/dx in place of 


tHe) — ii 
£2) = jin A 
dy d (dy 
7 ” i ees WS 
The symbol y" or ie is used for is (2). 
7 du. d (du 
and i or Ga is used for ai (#)- 


‘Third and higher derivatives are defined and written analogously. The nth derivative 
of f is denoted by f'"™. 


6.2 Rules of differentiation 


Sum rule: if wu and v are functions of x, then 
d du , dv 
a +v)= aie 
Constant multiplier rule: if k is a constant and u is a function of x, then 
(eu) = aoe. 
Product and quotient rules: if u and v are functions of x, then 
d du dv 
a = a” + uae 


e()=(e-*Z) /* 
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Chain rule or ‘function of a function’ rule: if f and g are two functions, then the 
derivative of their composition is given by the function of a function rule: 


d , 

ef (olz))1 = F'(G())9'(z)- 
A more common way to write this is the chain rule: 

dy _ dy du 

dx dudr’ 
where (in this case) u stands for g(x) and y for f(u). 
Implicit differentiation: given an equation connecting two variables, say 2 and y, 
we calculate dy/dx by differentiating both sides with respect to x and then solving the 
resulting equation algebraically for dy/dx (instead of solving for y before 
differentiating). 


6.3 Standard derivatives 


In each of the following cases the domain is the largest set of real numbers x for which 
the function is defined. 


Function Derivative 


ot azr®*-! (a any number) 
log, 1/x 

log.(—2) 1/x 

e e 

sine cos 

cos x sing 

tana sec? x 

cota —cosec? x 

secr sec rtanx 

cosecr —cosecrcotr 


The following derivatives are given in a form that is useful for integration. The 
constant a is assumed positive, and in the derivatives of arcsec and arccosec it is 
assumed that « > a (rather than x < —a). 


Function Derivative 


g 
3 
& 
5 
— 
ee] 


BIH BIR SIR SIH SIN 


Vw 


ae ee 


5 
i 
a8 
a 


6.4 Taylor series 


About x = 0: 
i og heey 
Mee ee Tete ak 
ee 
cose = 1 oe + ay ap tt 
2 
2 aoe 
faites ta tok 
2 
log.(1+2) = 2-545 —... (-l<2<1), 


=i 
(a+2)" =a" +ra"!e + ne Dgr-23 +: (-a<z<a) 


(where r is any real number and a is any positive number). 
In general, about 2 =0 
Sle) = $(0) + 2f'(0) + 5.27F"(0) + 225 F"(0) +>, 
and about ¢ =a 
f(x) = f(a) + (@ - a) f'(a) + Fle —a)?f"(a)+---. 
The nth-order Taylor approximation to f(x) about x =a is 
Sle) = fla) + (ea) f(a) + (0 ~a)?f"(a) +--+ A(x a)" fa). 
In particular, n = 1 gives the tangent approximation 
c J (x) = f(a) + (x — a) f(a) 
and n = 2 gives the quadratic Taylor approximation 
S(x) = f(a) + (@ — a) f'(a) + 3(@ — a)? f"(a). 


These approximations are good when - is close to a. 
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7 Integration 


7.1 Standard integrals 


The following integrals are obtained by reversing the entries in the tables of standard 
derivatives. In each case n stands for a non-negative integer, a for a non-zero real 
constant and a for a real constant. 


Function Integral Domain 
nyt 

a” ate All 

n+l 

goth 
2 (a¢-1) —S+e r>0 

(and x < 0 as well if a is an integer) 
x 
= log. 2 +C z>0 
2 
: bes(-2)40 =<0 
e tet +0 Alle 
1 

sinar Pt elaed +C Alle 

1 
cos ar an ar+C Alle 
sec? ax dtanaz + —r<ar<% 

a 2 2 
cosec? ax —* cotaz+C O<ar<nr 

1 T ™ 
secartanar —secar+C -~<ar<5 

a 2 2 


1 
cosecarcotaz ——cosecar+C O<ar<r 
a 


med 
arcsin— + C -a<r<a 
a 


Leese Alle 
a - 


1 r 
—aresee—+C xc>a 
a a 
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The following integrals may also come in useful. As before, the constant a must be 
non-zero, but b,a, 3,7 can take any real values (except where stated otherwise). 


Function Integral Domain 
1 i 
= 7 lose(ax +b) +E ar+b>0 
1 » 
as 7 eee(—ax — b) +E ar+b<0 
log. ar w(log, (ax) — 1) +C ar>0 
1 us 
tanar ceeded a <ar< 3 
cotar 7 08e(sin ax) +C O<ar<7 
secar * tog, (sec az + tan ar) +C -f<ar<t 
* tog, (cosec az — cot az) +C O<ar<n 
1 a+xr 
jae (+2) +0 |z| < Jal 
1 Z-a 
plvee(2=*) +0 z> {al 
1 —r+a 
je (=*2) +c z<-—la| 
log, (2 + Vz? + a2) +C All x 
log, (a + Vx? — a?) +C 2 >|a| 
log.(—a — Vz? — a?) +C xz <-t|al 


9 cer va 
(a >0) 
1 1 az —p 
= ‘COS +C 
V-aa? + 20x +7 va VF +ay 


(a > 0,6? + ay > 0) 


Vera 


z= a 


Va—a? 


2 
Flee e+ Va? + a2) + sve ra +C 
2 
Sloe, (2 — vz? —a?) + ave =a+C 
2 
F lee(—2 + Vr — a?) + sve =a#+C 


3 
P 


tog, (= + ai +3 +7) +C 


F wresin = + 5 Va + 


If ax? + 282 + 

has real roots, 

x > the larger root; if 
not, x is unrestricted. 
x between the roots 
of — ax? + 282+ 


Alla 
x > |al 
x <-—|la| 


2] < Jal 
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For the following integrals, the domain comprises all real values of x. 


Function Integral 


e* sin Br Sip lesin Be — Beos Be) +C 
a 


one 
€** cos Ba Sap pla cosie + Gsin Bx) +C 
a 


A x Dace 
rsinar ——cosar + — sinar+C 
a a 
z.. 1 
© cos ax ~sinaz + — cosar+C 
a a 
ate ae. 2 ar. 
wv sinax ——+- } cosar + —sinar+C 
aa a 
2 
@ eh ee 2a 
z* cosar (= - 5) sinar +B oosar +c 
a a a 
2 a n = 
x” sinax ——cosar + 2 fe ‘cosaxdx + C 
a 
ey (ile n seat 
2" cosaxr GZiner— 2 fs sina dr +C 
2 1 a 
sin* ax se — = sin2ar+C 
2 da 
r oe 
cos” aa st +> sin2ar+C 
2° 4a 


a 1 
sinarcosar —7- cos 2ar+C 


In the following three cases we require a # § and a # —3. 


Function Integral 


sin(a — f)r Es sin(a + 3)x 
2(a — 9) 2(a + 3) 
sin(a—B)r | sin(a + 3)x 

cos aw cos Ga We—8) + “2a+8) +C 
Gj _ cos(a —A)xr i cos(a + 3)xr 
sin ax cos Br We-6) ~ 2a+8) +C 


sin az sin Ox +C 


7.2 Rules of integration 


Examples 
1. The fundamental theorem of calculus 
[ fords = Fe) +e [rde= 4240, 
* é dF (xr) d 
if and only if f(z) = re r= Fra I} 


2. Linearity rule 
J tose) + da(e)) ax fi 
=a f seas +s f a(2)der 
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3. Substitution rules 


For integrals of the form To find J sitoz +6)dr =I, let 
t= [ase 42 dz, u=ar +6; then adz = du. So 
1 
l=- | si b)adx 
substitute u= f(x) and use Mae = di to ars J sn(ox+ a 
obtain 


= 2 fsinudu=—Heosu+c 
T= [ o(uydu. aJ a 
= 3 alan +b) +c. 


Otherwise, find a suitable relation between u To find / fesard = I, let 


eae alts is0 x =tanu; then dz = sec? udu. So 


jr 
dz =, a a sec? udu 
2* J + tan? us? 
sec? udu 


secSu 
= f cosudu =sinu +c 


= sin(arctan r) + C. 
4. Integration by parts 


Let u = u(a) and v = v(x); then To find J rcoszde =I1y, letu=a 
[vGar=w- f and v = sinz; then 
dr 


is = [ eZsinz) az 
=arsinz— / snz x ldr 


=arsinz+cosr+C. 


7.3 Integration methods 


Rational functions 
To integrate a rational function, that is, one of the form P(r)/Q(x) where P(x) and 
Q(x) are polynomials, we first express the function as a sum of partial fractions (see 
Subsection 3.5). Then integrate each partial fraction as follows: 
dx 
ar+b 


a standard form; 
/ mF substitute u = ax + b; 


ott de with 6? < dac, 
eee) substitute 2ar + b = Vac — 6 tanu. 


Trrational functions 
Integrals containing the expressions a? — x?, Vx? — a?, Vx? +a? can often be 
simplified by the following substitutions: 
for Va?—2? use x =asinu; 
for Vx —a@ use x =asecu; 
for Vrt+a2 use x©=atanu. 
More generally, for y/ax? + Bx +7: 
ifa>OandB?<4ay use 2or+B=V4ay—# tanu; 
ifa>Oand 8? >4ay use 2ar+B= VF —4aysecu; 
ifa<Oand#>4oy use 2ar+3= fF —4aysinu. 
(In the remaining case the square root is always imaginary.) 


Exponential and trigonometric functions 


For if P(x) {im ar \ dx, where P(x) is a polynomial, integrate by parts (several times 
cosar 


if necessary) with u = P(2), = =e” or sinaz or cosazx. See also the ‘standard 


integrals’ table for if x" sinardx and / 2” cosar dx. 


ax fsin Bar ‘ * + 
For fe ee {ene} an, see the ‘standard integrals’ table. 


For i sin" ax cos" ax dx, 


ifnisodd use u=sinaz; 
ifmisodd use u=cosar; 
if m,n are both even, then 
ifn<0 use u=tanaz, 
ifm<0 use w=cotar, 
if m,n both >0 use sin?ar = $(1 —cos2az), cos*ax = }(1 + cos 2az). 


For other integrals involving trigonometric functions only, it may be helpful to express 
them in terms of sinax and cosaz and then use the substitution t = tan }az, according 
to which 

-? 2dt 
pes ad: 


4 2t 
sinat = cosaz = 77 dee I 


1+t 
2 


1-# 


For integrals involving e* only, the substitution u = e** may be useful. 


8 Unit outlines 


Unit 1 Recurrence relations 


Section 1: Introduction to recurrence relations 


1. An equation relating each term in a sequence to one 
or more previous terms, which can be used to generate 
the whole sequence once the first few terms are known, is 
called a recurrence relation. 


2. A recurrence relation is of kth order if the difference 
between the highest and lowest subscripts in it is k. To 
generate a sequence using such a recurrence relation we 
need in addition the values of the first k members of the 
sequence; any set of conditions from which they can be 
deduced are called the initial conditions. 


3. A recurrence relation is linear if it can be written as 
Ung) = pty + betlnnt + Crtp—a boo + Pry 
where p, and the coefficients a,, b,, cr, ... may depend 


on r, but do not depend on any of the us. Otherwise it is 
said to be non-linear, 


4. A linear recurrence relation has constant 
coefficients if none of the coefficients ay, by, cr, . 
depend on r. (But p, may depend on r.) 


5. A linear recurrence relation is homogeneous if 
pr =0 for all values of r. Otherwise it is said to be 
non-homogeneous (or inhomogeneous). 


6, The general solution of a recurrence relation is a 
formula containing one or more arbitrary constants, from 
which all the particular solutions can be obtained by 
giving the arbitrary constants particular values which 
may be determined from the initial conditions. 


7. The general solution of the linear 
constant-coefficient first-order recurrence relation 


Upp = ay +P 


Ba” 
un = 
A+np 


where A and B are arbitrary constants. 


is 


a=l, 


Section 2: Linear second-order recurrence 

relations 

To obtain the general solution of the linear homogeneous 

constant-coefficient second-order recurrence relation 
Ure = aur + buy, 

we first consider the auxiliary equation 
zr =ar+b. 

There are three cases to consider_ 

(i) Ifa? +46 > 0, the auxiliary equation has two 
distinct real solutions, and 1, and the general 
solution for the recurrence relation is 


uy = AN" + Bu". 
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(ii) If a? + 46 = 0, the auxiliary equation has a repeated 
real solution \, and the general solution is 


Un = (A+ Bn)d". 


(iii) If a* + 4b < 0, the auxiliary equation has no real 
solutions. (See Unit 5, Section 5.) 


Section 3: Numerical difficulties 


1. If Z is an approximation to some quantity whose true 
value is x, we call  — x the absolute error in z, and we 
call (% — x)/x the relative error in x. 


2, Consider any problem in which we are given numerical 
data and we want to calculate some result from them. 
The problem is said to be absolutely ill-conditioned 
(with respect to small changes in the data) if a small 
absolute error in the data gives rise to a significantly 
larger absolute error in the result. It is said to be 
relatively ill-conditioned if a small relative error in 
the data gives rise to a significantly larger relative error 
in the result. It is said to be absolutely 
well-conditioned if the absolute error in the result is 
smaller than the absolute error in the data, and 
relatively well-conditioned if the relative error in the 
result is smaller than the relative error in the data. 


3. The seale factor of a calculation is the factor by 
which a small error in the data is multiplied to give the 
corresponding error in the result. If the magnitude of the 
scale factor is significantly greater than 1, then the 
problem is absolutely ill-conditioned with respect to small 
changes in the data; if the magnitude of the scale factor 
is less than 1, then the problem is absolutely 
well-conditioned. 


4. Consider the problem of calculating un from uo, using 
the recurrence relation 

Urs1 = aur +p. 
The scale factor for this problem is a”. If the magnitude. 
of a” is less than 1, the problem is both absolutely and 
relatively well-conditioned with respect to small changes 
in uo. If the magnitude of a” is significantly larger than 
1, then the problem is absolutely ill-conditioned. For any 
a such that ja! > 1, the problem is therefore absolutely 
ill-conditioned for sufficiently large n. It is relatively 
ill-conditioned if, in addition, the condition 


Up = 


is satisfied. 


5. If |a-| > 1 for all r, then the problem of calculating up 
from uo using the more general recurrence relation 

Upp1 = Ary + Pr 
is absolutely ill-conditioned (with respect to small 
changes in uo) for sufficiently large n. It is relatively 
ill-conditioned if, in addition, the initial value uo is close 
to a value for which u,,, the true solution of the 
recurrence relation, remains small for large values of n. 


Section 4: Living with ill-conditioning 


1. Mortgage problems, relating to a situation where £X 
are borrowed at J% per annum interest and repaid by 
monthly payments of £M, can be modelled using the 
recurrence relation 


urg1 = (1 +1/100)uy — 12M, 
where u, is the amount owed at the end of the rth year. 
The initial condition is uo = X. 
Mortgage problems are often highly sensitive to small 
changes in the data, particularly the values of J and M. 


2. Sometimes ill-conditioning can be cured by a 
reformulation of the problem. For example, in the 
recurrence relation 

Urt1 = Orly + Pry 
with |a,| > 1 for all values of r, the solution which 
remains small cannot be reliably calculated by forward 
recurrence, i.e. starting with up and calculating u:, 
uz, ... in succession, because this problem is relatively 
ill-conditioned, The reformulation which cures the 
ill-conditioning is to use backward recurrence instead. 
This means rewriting the recurrence relation with uy as 
the subject of the formula, 

t= SRS, 

ce 

and calculating values of uy from it working down from a 
large value of r, say R. Even though we know very little 
about the value of up, it is easy to obtain accurate 
results because the new problem is very well-conditioned. 


Unit 2 Differential equations I 
Section 1: Direction fields 


1. A first-order differential equation is an equation 
connecting the derivative dy/dx of some unknown 
function to the variables x and y which are related by 
that function. By solving the equation algebraically for 
dy/dx we can bring it to the standard form 


dy _ 
ae me) Q) 
where m(z,y) stands for some formula which may 


contain either or both of the variables z and y. 


2. By a solution of a differential equation we mean a 
function—that is, a rule or formula giving y in terms 
of z—which satisfies the differential equation. 


3. A direction field is a rule or formula which specifies 
a slope at each point in the (.r, y)-plane (or some part of 
it). To each direction field there corresponds a differential 
equation, constructed by taking m(:,y) in (1) to be the 
formula giving the slope of the direction field at the point 
(x,y). If y= f(z) is the equation of a trajectory of the 
direction field—that is, a smooth curve at every point of 
which the slope of the curve matches that of the direction 
field—then y = f(z) is also a solution of the 
corresponding differential equation. The direction field 
has an infinite family of trajectories, each one 
corresponding to a different solution of the differential 
equation, 


4. For some differential equations, an exact formula can 
be found for the solutions. This formula, called the 
general solution, contains an arbitrary constant of 
integration, often denoted by C. Different numerical 
values for C give different solutions to the differential 
equation (usually called, in this context, particular 
solutions to distinguish them from the general solution) 
and therefore label different trajectories of the direction 
field. 


5. Whether or not it is possible to find a formula for the 
solution of a differential equation, the direction field is a 
useful source of information about its solution; in 
particular it can be used to sketch solution curves 
(trajectories) and so obtain an idea of their qualitative 
behaviour. 


Section 2: A numerical method 
Given a differential equation 

dy 

<= m(x,y) 


dr 
and the condition that y(z) = yo, where xo and yo are 
specified, then Euler’s method can be used to find y(b), 
where 6 is a given value of x, as follows. 
(i) Choose h = (b~ 9)/n, where n is the number of 
steps to be used, 
(ii) For r= 0,1, n—1, use the recurrence relation 
Your = ¥e + hm(ay, Ye), 
where 2, = to + rh, so that x, = b, and Yo = yo. 


(iii) Then Yp is an approximation to y(b). 


‘The error in the approximation to y(b) is roughly 
proportional to h, and therefore inversely proportional to 
the cost (i.e. number of steps) of the calculation, 


Sections 3 and 4: Exact methods 


1, To look for an exact solution to a given differential 
equation, first bring it to the standard form 


® = (x,y). 

Then 

(i) if m(x,y) has the form 
m(x,y) = f(x)g(y) 


(where either or both of f(z) and g(y) may be 
constants) use separation of variables (see below); 


(ii) if m(,y) has the linear form 
m(x,y) = k(x) + U(r)y 
(again either or both of k(x) and U(r) may be 
constants) use the integrating factor method (see 
below). 
2. Separation of variables (Procedure 3.3) 
To solve the differential equation 


dj 
= s(x)a(u) (alu) #0), 
divide both sides by g(y) and integrate, to obtain 


1 
[ape [soure. 


Perform the two integrations and make y the subject of 
the resulting equation. 


3. Integrating factor method (Procedure 4.2) 
To solve the linear equation 


2 = k(z) + Ux)y, 


subtract [(x)y from both sides and then multiply both 
sides by the integrating factor 


P(x) = exp [[cuenae] : 


The resulting equation has the form 


pte) de ap 2), =4(x) 


and integrates directly to give 


r(x)y = [oo dr+C, 
which is easily solved for y to give the general solution. 


4, Given the general solution of a differential equation 
involving an arbitrary constant, say C, we can determine 
the particular solution satisfying a given further 
condition by substituting the general solution into the 
further condition. This gives an algebraic equation which 
we can solve for C. Inserting this value of C into the 
general solution of the differential equation gives the 
required particular solution. 


5. If you need help with the integrations, consult the 
Integration section of this Handbook. 


Unit 3 Animal populations: their 
growth and exploitation 


Section 1: How to model populations 


Mathematical modelling is the use of mathematics in 
solving real-world problems. Empiricism means basing 
one’s model on experiment, without any theory; in the 
study of population changes it means fitting a simple 
mathematical relation to the data without seeking to 
explain why it should have that particular form. If the 
empirical relation is then applied outside the range 
covered by the original data, we are extrapolating, a 
notoriously inaccurate procedure. 
A useful tool for setting up models which do have some 
theoretical basis is the input-output principle: 
increase in population (i.e. accumulation) 
= input — output. 
The exponential model is obtained by assuming that 
(i) the population P(t) can vary continuously, 
(ii) migration is negligible and (iii) the numbers of births 
and of deaths in any short time interval of duration A are 
respectively BP(t)h and CP(t)h, where B and C, the 
proportionate birth rate and proportionate death 
rate, are taken to be constants. The resulting equation is 
dP 
a = (B-O)P. QQ) 
The long-term predictions of this model are 
(i) unbounded growth if B > C, (ii) decay to zero if 
B <C and (iii) constant population if B = C. Case (i) in 
particular is clearly unrealistic, but for some populations 
these predictions fit the real data reasonably well over 
limited time intervals. 


Section 2: The logistic model 


In the logistic model we assume that B and C, instead 
of being constants, depend on the total population in 
such a way that B — C, the proportionate growth 
rate, decreases linearly with the population P. The 
resulting equation is the logistic equation 

@~*(-77)* ) 
where a and M are constants. Built into this model are 
the assumptions that 


(i) there is an equilibrium population level, M; 


(ii) the proportionate growth rate depends only on the 
instantaneous value of P; 


(iii) the proportionate growth rate decreases linearly with 
P, being given by a(1 — P/M). 
The solution of the logistic equation is 


M 
ETS To (M/P(0) — tew*** 


and the graph of this solution (for 0 < P(0) < M) is 
called a logistic curve. 


t 


‘Two ways of comparing this model with observed values 
for P(t) are: 


(i) estimate M from the property that P(t) is close to 
M for large t, then plot log.(M/P(t) — 1) against t; 


(ii) if data are available at equally-spaced instants 
ty, ta,.-+, plot (Pay1 — Pa)/Pn against Pair. 


In each case the logistic model predicts that the graph 
will be a straight line. 

Though far from being a ‘law of nature’, the logistic 
model! agrees quite well with data for some populations, 
particularly under laboratory conditions; but for wild 
populations the agreement is more likely to be qualitative 
than quantitative. In cases where the model does not fit 
the data well, it is usually possible to see that one or 
other of the assumptions (i) to (iii) above is not realistic. 


Section 3: Fishing 
The following procedure, known as the Schaefer model, 


can be used to estimate optimum fishing effort and 
maximum sustainable yield for a fishery. 


(i) Plot yield per unit effort (a measure of population) 
for each year against the mean fishing effort over the 
preceding few years. 

(ii) Fit a straight line to the points plotted in (i) and 
assume this to give the relation between the 
long-term sustainable yield per unit effort and 
sustained fishing effort. 


(iii) Deduce the relation between sustainable yield and 
sustained fishing effort, and hence find the sustained 
fishing effort that gives the maximum sustainable 
yield. 
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The Schaefer model assumes that 

(i) the modified logistic model (3) below applies; 

(ii) the historical variations of fishing effort were slow 
enough for the fish population in each year to be 
close to the equilibrium level appropriate to the 
mean level of fishing effort exerted over the previous 
few years; 

(iii) catch is proportional to the product of fish 
population and fishing effort. 


Section 4: Models of exploitation 


‘We can modify the logistic model to include the effect of 
exploitation by human beings. The new equation is 
dP jel 
Peaa(i-7)P-c, (3) 
where C is the rate at which members of the population 
are caught, i.e. it is the yield. Two exploitation 
strategies are considered. 
(i) Constant catch, i.e. C = constant. Here the 
conclusions are 
(a) if C exceeds the maximum sustainable yield 
aM/4, the population will dwindle to zero; 
even if C <aM/4, the population will dwindle 
to zero if P(0) is too small; 
(c) if C <aM/4 and P(0) is large enough, then for 
long times P(t) will approach an equilibrium 
value smaller than M. 


(b 


(ii) Constant effort, i.e. C proportional to P, say 
C=AfP, where f (the effort) and \ are constants, 
Here the conclusions are 
(a) if Af > a, then the population will dwindle to 
zero; 

(b) if Af <a, then the population is viable whatever 
the value of P(0); 

(c) in the long run, maximum yield is obtained with 
a value of Af = 4a. 


A more general model, in which (3) is replaced by 

dP 

lt Act fe 
with f(P) a suitable function of P, gives qualitatively 
similar conclusions, The qualitative features of the above 
conclusions are therefore more reliable than the 
quantitative ones. 
An important difference between the two strategies is 
that if we aim for maximum yield, the constant catch 
strategy is unstable—a small variation in the conditions 
can cause a large change in the result, and can in this case 
cause the population to dwindle to zero—whereas the 
constant effort strategy is stable—a small variation in 
the conditions produces only a small change in the result. 


Unit 4 Newtonian mechanics in one 
dimension 

Section 1: The description of motion 

1. A particle is a material object whose size and 


internal structure are negligible, so that at any given time 
it is located at a single point. 


2. The motion of a particle along a straight line is 
described by a position function z(t). For each instant 
t in time, this function specifies the position x(t) 
occupied by the particle. 


3. The velocity v of the particle is the derivative of the 
position function with respect to time, that is, 
va 
dt” 
4. The acceleration a of the particle is the derivative of 
the velocity function with respect to time. In symbols, 
dv dx du 


ze 


ade de 


5. In the special case of uniform motion, by which is 
meant motion with constant velocity uo, we have 


x2 =zo + vot, 
v=o, 
a=0, 


where zo is the position at time t = 0. 


6. In the special case of motion with constant 
acceleration ap, we have 

x=z0 + vot + Sot”, 

v = vo + aot, 

v? = vg + 2ao(x — x0), 
where xo and vp are the position and velocity of the 
particle at the instant t = 0. 


Section 2: Newton’s laws of motion 


1. Newton's first law states that, in the absence of a 
force, a particle either stays permanently at rest or moves 
at constant velocity. 


2. The force acting on a particle at a given instant is 
represented by a single number, F’, called the force’s 
z-component. The force has magnitude |F|, and is in 
the direction of increasing x if F > 0 or in the direction 
of decreasing x if F <0. Often ‘the force with 
z-component F” is shortened to ‘the force F’. 


3. Newton's second law states that if a particle of 
mass m experiences a net force with z-component F, 
then its acceleration a is given by F = ma. 


4. If more than one force acts on a particle, then the net 
force or total force F to be used in Newton's second 
law is the algebraic sum of the individual forces, that is, 
the sum of their z-components. 


Section 3: First model for the motion of a 
falling object 

In this section Newton’s second law is used to model the 
motion of an object falling under the force of gravity, 
where the effects of air resistance are ignored. 


1, Newton's second law F = ma may in principle be 
integrated by putting 


(i) 


2, Any object near the Earth's surface is pulled 
downwards by the force of gravity. The magnitude of 
this force is proportional to the object’s inertial mass m, 
that is, 

force of gravity = mg downwards, 


where the constant of proportionality g is approximately 
9.81ms~*. 


Section 4: Second model for the motion of a 
falling object 


In this section Newton’s second law is used to model the 
motion of an object falling under the forces of gravity and 
air resistance. 


1. The motion of any object through the Earth's 
atmosphere is opposed by the force of air resistance, 
whose magnitude R depends on the object's speed |v}, 
shape and size. For small speeds it is appropriate to use 
the linear approximation 
R= kilo), 
whereas for larger speeds the quadratic approximation 
R~ kalv|* 
is appropriate. Here k; and k2 are constants which 
depend on the object's shape and size. 


2. The constant speed |u| at which an object can fall is 
called the terminal speed. Using Newton's first law, 
this occurs when R = mg. 


Unit 5 Complex numbers 


Section 1: Introducing complex numbers 


The natural numbers are 1, 2, 3, ...; the integers are 
these together with their negatives and zero; rational 
numbers are numbers of the form m/n, where m, n are 
integers (n #0); real numbers can be thought of as 
possible lengths of line segments and are often 
represented geometrically as points on a number line. 
An irrational number is one that is real but not 
rational, 
A complex number can be defined as a pair of real 
numbers (x,y) obeying the following laws of addition and 
multiplication: 

(1, yx) + (wa, y2) = (a1 +22,y1 + y2), 

(1,91) x (#2, y2) = (122 — yaya, tiys + yi22). 
Like the rational and real number systems, the complex 
number system is arithmetically complete in the sense 
that the sum, difference, product and quotient of any two 
numbers in the system (with the single exception of 
division by zero) is again a number in the system. 
If 2 is a complex number (z,y), then the real part of = is 
Re: =z and the imaginary part is Im z= y. Two 
complex numbers are equal if and only if their real parts 
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are equal and their imaginary parts are equal. A complex 
number is (purely) real if its imaginary part is zero, and 
(purely) imaginary if its real part is zero. 
The following notation, called the Cartesian form, is 
customary for complex numbers: 

(x,y) is written 2 + iy, 

(x,0) is written x, 

(0, y) is written iy or yi, 

(0,1) is written i, 

(0,0) is written 0. 
Calculations with complex numbers are then easily 
carried out by following the usual rules of algebra with 
the additional rule i? = —1. A list of these rules is given 
in Subsection 3.6 of this Handbook. 
The complex conjugate of a complex number 
z= 2+ iy is the complex number = = x — iy, which has 
the property that 23 is real and equal to z? + y?. To 
calculate the quotient of two complex numbers, multiply 


numerator and denominator by the complex conjugate of 
the denominator. 


Section 2: Geometrical representation 


The Argand diagram is a representation of the complex 
number system, in which the complex number x + iy is 
represented by the point in the plane with Cartesian 
coordinates (x,y). Addition of complex numbers can then 
be represented by a parallelogram construction, and 
complex conjugation corresponds to reflection in the 
a-axis. 


Section 3: Polar form 
If the point representing a complex number z has polar 
coordinates [r,6] (see Subsection 4.4 of this Handbook), 
we can write = in the polar form 

= =r(cos@ + isin@). 
‘The modulus of z is || = r, and 6 is called an 
argument of =. If @ lies in the principal value range 
—7 <@< 7x, we call @ the argument of = and write 
8 = Arg z. In terms of x and y, the real and imaginary 
parts of z, we have 


l2]= V2? +y", 
ify >0, 
Arg: = { *eeos(z/|2)). ify 20 
—arccos(a/|2|) ify <0. 
When two complex numbers are multiplied, the moduli 
are multiplied and their arguments are added. In 
particular, de Moivre’s Theorem states that, for any 
integer n, 
(cos @ + isin @)" = cos nd + isin nd. 


Section 4: Euler’s formula 
Euler’s formula 
ef” = cosO +isind 


enables us to write a complex number z in exponential 
form as 


e 


z 


where r= |2| and @ = Argz, and to define the 
exponential function of z as 


=e*(cosy+isiny), 
where x = Rez and y = Imz. 


expz=e=ett 


Euler's formula also gives the results 
1 SK sly er ae 
cos = 5(e" +e), sind = Ae = oa) 
which can be used to derive trigonometric identities. 
A sinusoidal function of time t is one having the form 
f(t) = Acos(wt + 4), 
where the constant A is called the amplitude, w the 
angular frequency and ¢ the phase. The function 
repeats after a time T’ = 27/w, called the period. The 
complex number Ae" is called the phasor of the 
sinusoidal function. Sinusoidal functions with the same w 
can be added by adding the corresponding phasors. 


Section 5: Complex roots 


If a = Re®® is a complex number and n is a positive 
integer, the nth roots of a, written a'/", are the 
solutions of the equation 2" =a, There are precisely n of 
them, given by 

all” = RY" exp (iP 47") (k= 0,1,....0=1)- 
The fundamental theorem of algebra states that any 
polynomial of degree n, 

P(z) =anz" +@n—12" + +++ +412 +40 
(with an # 0), can be factorized into n factors of degree 1, 

an(z — 21)(2 — 22)... (2 — =n); 
where z1,..-,2» are complex numbers. These numbers 
are called the roots of the polynomial and are the 
solutions of the equation P(z) = 0. If the coefficients 


Qn, Qn—1,.++,@0 are all real, then the non-real roots (if 
any) occur in complex conjugate pairs. 


For second-order recurrence relations of the form 


n=1 


Ur¢2 = aury1 + bur, 


if the auxiliary equation z* = az + 6 has non-real 
solutions re*"®, the general solution is 


Un = r"(Acosnd + Bsin nd). 


Unit 6 Differential equations II 


Introductory definitions 


A second-order differential equation is an equation 
relating the second derivative of an unknown function to 
one or more of the following: known functions, the 
unknown function, and the first derivative of the 
unknown function. The general solution of such an 
equation normally involves two arbitrary constants. If a 
second-order differential equation cari be written in the 
form 
2 

pla) £4 + a(e) $2 + ray = 2), a) 
where p, q, 7 and f are given functions (and p() is not 
the zero function) then it is linear. If in addition f is the 
zero function, then the equation is homogeneous; 
otherwise it is inhomogeneous (or 
non-homogeneous). If a linear second-order differential 
equation has the form 


ot 4 6 + ey= slo), (2) 


where a, 6, c are constants with a # 0, then it is said to 
have constant coefficients. 


Section 1: Homogeneous equations 
To solve a homogeneous second-order linear equation 
with constant coefficients, i.e. one of the form 


age tba tey=0, (3) 
apply the following procedure. 
(i) Write down the auxiliary equation 
a’ +bA+e=0, (4) 
and solve it algebraically for A. 
(ii) (a) If (4) has two real distinct solutions A, and A2, 
then the general solution of (3) is 
y(z) = Ae™* + Be?*, 
where A and B are arbitrary constants. 
(b) 


If (4) has only one real repeated solution A; (i.e. 
a\* +6\ +c is a perfect square), then the 
general solution of (3) is 


y(z) = (A+ Baje**. 
If (4) has complex conjugate solutions 
Ay =a +i and Ag = a — iG, then the general 
solution of (3) is 


y(z) = e**(Acos Gr + Bsin fx). 
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If the differential equation is homogeneous but does not 
necessarily have constant coefficients, the following 
theorems tell us something about the solutions. 

Theorem 1.1: If the functions y; and yz are solutions of 
a homogeneous linear differential equation, then any 
linear combination of y: and yz is also a solution. 
Theorem 1.2: If, in addition, the differential equation is 
of second order, i.e, of the form 


2 
ple) saz) +r(z)y=0 

with the coefficient p(x) never taking the value zero, and 

if the solutions y; and y2 are linearly independent, then 

the general solution of the differential equation is 

Ay: + Bya, where A and B are arbitrary constants. 

By a linear combination of two functions y; and y2 

with the same domain we mean a function of the form 

ciy: + ¢2ya, where c; and c2 are constants. The functions 

yj: and yp are said to be linearly independent if neither 

is a constant multiple of the other. 


Section 2: Inhomogeneous equations 


For the general inhomogeneous Equation (1), we define 
the associated homogeneous equation to be 


wa) 4 + a(e) 2 + r(2)y=0. (6) 
The general solution of (5) is called the complementary 


function of (1), and is denoted here by ye. The general 
solution of (1) is (Theorem 2.1) 


Y= Vet Up 
where yp is any particular solution of (1). 
For certain constant-coefficient inhomogeneous equations, 
particular solutions can usually be found by substituting 
the function for yp suggested below into the differential 
equation and choosing the constants m, n or 2 (complex) 
so that the equation is satisfied. 


f(z) Try for yp(z) 
constant constant 

ka +l mr+n 

ke me” 
keoswa + lsinwz | mcoswx +nsinwa 
or Re(k —il)e““* | Re ze‘** 


In the above, k, |, a and w denote given constants. 
Exceptionally, the differential equation may have no 
solution of the suggested form; in that case, multiply the 
suggested function by x and try again. 

‘The principle of superposition: If u(r) is a solution of 
Equation (1), and if (2) is a solution of the equation 
obtained by replacing f(x) in (1) by g(x), then 

au(x) + bu() (with a and 6 constants) is a solution of the 
equation obtained by replacing f(x) in (1) by 

af(x) +bg(zx). 


Section 3: Initial and boundary conditions 


In an initial condition problem for Equation (1), the 
values of y and its derivative are both specified for the 
same value of 2. Such a problem always has a unique 
solution (so long as p(2) never takes the value zero). 


In a boundary condition problem, one condition is 
placed on the value of y (or its derivative) at each of two 
different values of x. Such a problem may have a unique 
solution, or no solution, or an infinite number of 
solutions. 


Section 4: Nature of solution functions 


For the homogeneous Equation (3), with a and c positive, 
the solution y() may depend on « in various ways, 
depending on the value of b: 


b< -2/ae (for large x) increasing exponential; 
—2/ae <6 <0 oscillation with increasing amplitude; 
b=0 sinusoidal oscillation; 

0<b<2Vae oscillation with decreasing amplitude; 
2Vac <b (for large x) decreasing exponential. 

For the inhomogeneous equation 
2 
att + ot +cy = kcoswr + Isinwz, 


with a, b, c, k, |, w constants and a, b, ¢ all positive, the 
complementary function decays to zero as x increases, 
and is called a transient. For sufficiently large x, any 
solution of the equation is close to the unique purely 
sinusoidal particular solution, which is called the 
steady-state solution. 


Section 5: A numerical method 


Any second-order differential equation 
@y dy 
ftom (anit 

can be written as two simultaneous first-order equations 
dy 
dx 
© = m(x,y.2)- 

Euler’s method approximates these by the recurrence 

relations 
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You1 = ¥, +hZ,, 
Zr41 = Z, +hm(z,r,¥;,2Z,), 


where h is the step size, rz, = Zo + hr and Y,, Z, are the 
approximations to y(z,) and 2(z-). 


Unit 7 Oscillations and energy 


Section 1: A home-made oscillating system 
Section 2: The perfect spring 


In this section a perfect spring is used to model vibrating 
systems. 


1, A perfect spring is characterized by two positive 
constants: its natural length and its stiffness. It obeys 
Hooke’s Law, exerting the following force on any object 
attached to either of its ends: 


(i) when the spring is eztended, the force is directed 
towards the centre of the spring and has magnitude 
(called the tension) equal to stiffness x extension; 


(ii) when the spring is compressed, the force is directed 
away from the centre of the spring and has 
magnitude (called the thrust) equal to stiffness x 
compression. 


2. Ifa particle in one-dimensional motion is attached to 
perfect springs then the equation of motion of the 
particle can be derived by considering just one 
configuration in which each spring is designated as being 
either extended or compressed. 


3. The equation of simple harmonic motion is 
#+w*r=0, 
whose general solution can be written either as 
z= Beoswt + Csinwt 
or as 
«= Acos(wt + 6). 
The constant w is called the angular frequency and the 
constant @ the phase angle of the oscillations. The 
period of the oscillations (time for one complete cycle) is 
7 =2n/w. The frequency (number of cycles per second) 
is f = 1/r=w/2n. The amplitude of the oscillations is 
A. The constants A, ¢ and B, C are related via the 
equations 


A= VB?+C? 
and 


_f arccos(B/A) C <0, 
= —arccos(B/A) C>0. 


4. For a particle in simple harmonic motion, the 
equation = +w*z = 0 describes the motion when the 
origin of the x-axis is chosen to be at the equilibrium 
position of the particle. If another origin is chosen then 
the equation of motion becomes 


f+Wr=wx,, 
where x = 2, is the particle's equilibrium position. For a 
particle of mass m attached to one end of a perfect spring 
of stiffness k whose other end is fixed, the angular 
frequency w is given by w? = k/m for either vertical or 
horizontal motion. Here gravity is taken into account in 
the vertical case, but other forces are ignored. 
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Section 3: The conservation of energy 


1. The relationship between a force with x-component 
F(z) and the potential energy function U(z) is 


U(z)= ~ f Flea, 
or, equivalently, 
F@)== ee 
The point at which the potential energy function is zero 


is called the datum. This can be chosen to be any 
convenient point. 


2. In particular, the gravitational potential energy is 
U = mg x height, 
where the height is measured from the chosen datum. 


3. The potential energy of a perfect spring is 
U = } x stiffness x (extension)*. 
4. The law of conservation of mechanical energy 
states that if the total force acting on a particle depends 
only upon the particle’s position then the quantity 
dmv? + U(x) = EB 
remains constant throughout the particle's motion. The 
quantity 
T = hmv" 
is known as the kinetic energy of the particle, and the 
constant E as its total mechanical energy. 


5. The region of motion of a particle with total 

mechanical energy E satisfies the inequality 
B-U(z)>0. 

‘The end-points of this region, for which 
E-U(x)=0, 

are called the turning points of the motion. 


Section 4: The fabulous perfect spring 


For a system such as a pendulum or a toy car on a curved 

track, the law of conservation of energy has the form 
}m@’ + U(a) = B, 

where q is the displacement of the pendulum bob along 

its path or of the car measured along the track. For 

small oscillations about the lowest point of the path, 

q = qo, the potential energy function can be 

approximated by the first three terms of its Taylor series. 

This approximation leads to the prediction of simple 

harmonic motion of angular frequency w ~ \/U"(qo)/m 


and period r= 2ny/m/U"(qo)- 


Unit 8 Damped and forced vibrations 
Section 1: Damped vibrations 


1, Damping car be modelled by a damping force of 
magnitude r|é|, where r is a positive constant and z is 
the velocity of the particle, and direction opposed to the 
velocity of the particle. The SI units of the damping 
constant are Nm~'s (or kgs~'). 


2. The equation of motion of a particle of mass m 
moving under the influence of a perfect spring of stiffness 
k and linear damping with damping constant r is 

mi +rz+kz=0, 
where zx is the displacement of the particle from its 
equilibrium position. Such a system is called a damped 
harmonic oscillator. 


3. For weak damping (r? < 4mk) the general solution 
of the equation of motion is 


a(t) = Ae~* cos (Qt + 4), 
where 


eas oo ane 


All solutions of this form may be obtained with the 
arbitrary constants A and ¢ restricted so that A is 
non-negative and ¢ lies between —x and 7. 

This motion is oscillatory with period r= 2x/®. The 
amplitude of the oscillations decreases exponentially, so 
that 


x(t + 7)/2(t) = e~ 


pr 


4. For critical damping (r* = 4mk) the general 
solution is 

a(t) = Be" + Cte, 
where p = r/2m and B and C are arbitrary constants. 
This case gives the fastest return to the equilibrium 
position. 


5. For strong damping (r? > 4mk) the general solution 
is 


x(t) = Be?!" + Ce", 


where 
r-VP=amk |, tv amk 
a= am Pips mm 


Section 2: Forced vibrations 


1. The equation of motion of a damped harmonic 
oscillator subject to an applied sinusoidal force of 
amplitude P and angular frequency w can be written as 


mz +rz+kx = Pcosut. 
The steady-state forced vibrations of this system can 


be found by using the phasor method. 
2. This equation of motion can be rewritten in the form 
# + 2awot + wor = £ coswt, 


where the undamped angular frequency w is defined 
(as in Unit 7) by 


k 
w= —. 
m 
and the damping ratio a is defined by 
, 
~ 2Vmk 


3. The steady-state vibration of the forced and damped 
harmonic oscillator is given by 


x = Acos(wt +4), 


where 
eee. Ss 
© lak = wo) + da2akat 
and 


A 
ate eee 
Vu? — 0)? + 402 


4. For a < 1/2, the amplitude of the forced vibrations 
exhibits a maximum at a particular value of the forcing 
angular frequency w. ‘This effect is known as resonance. 


Section 3: The perfect dashpot 


1. A perfect dashpot represents a resistive force which 
is proportional to the relative velocity between two 
components of a mechanical system. 


2. The dashpot force has magnitude R = r{/|, where I is 
the length of the dashpot and r is a positive constant 
known as the dashpot constant. If the dashpot’s length 
is increasing, then the dashpot force is directed towards 
the centre of the dashpot; if the length is decreasing, then 
the force is directed away from the centre of the dashpot. 


3. The equation of motion of a particle attached to a 
perfect dashpot can be derived by assuming that the 
dashpot's length is increasing. This equation of motion 
applies also to the case where the dashpot is contracting. 


Unit 9 Simultaneous linear algebraic 
equations 


In this unit the methods and theory are described for the 
special case of 3 equations in 3 unknowns, but the 
descriptions can easily be extended to the case of n 
equations in n unknowns. 

A set of 3 simultaneous linear algebraic equations 
in the 3 unknowns 21, x2, 73 is one that can be written 
in the form 


grmthrtears=dh, Fy 


aga, +borot+ers=d2, Ey 
a3z, + bara + egrt3 = ds. EBs 


The array of numbers 
a bh 
az ba cy 
ay bs cy 


is called the matrix (of coefficients) of the set of 
equations, and the array 


a bh a | dy 

[ees “| 
ay bs ¢3 | ds 

is called the augmented matrix. 


Section 1: Gaussian elimination 


A systematic procedure for solving such sets of equations 
is Gaussian elimination. 


Stage 1: Elimination 


(i) Eliminate 2; from all but equation Ey, by 
subtracting suitable multiples of E; from each of 
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the equations below it. In the above example the 
multipliers used in calculating these multiples of E; 
are a2/a; and a3/a; for equations E> and Es 
respectively. 

{ii) Eliminate x2 from all but the first two of the new set 
of equations by subtracting a suitable multiple of the 
new second equation from the equation below it. 


(iii) The new set of equations now has the ‘upper 
triangular’ form 


a2, +br2+ crs = di, Ey 
b322 + ¢}23 = d3, Ex. 
rs=d3. Bay 


Stage 2: Back-suzbstitution 
(i) Solve Ey, for x3. 


(ii) Substitute this value of xy into B2,, and solve the 
resulting equation for x2. 


(iii) Similarly, substitute these values of x3 and x» into 
E\, and hence find x}. 


Section 2: Types of solution 


1. A set of linear simultaneous equations will have: 
(i) a unique solution, or 
(ii) no solution, or 


(iii) an infinite number of solutions. 


To find out which case applies, look at the augmented 
matrix of the set of equations produced by stage 1 of 
Gaussian elimination. [In the table below, a matrix entry 
which must be non-zero is denoted by v’, while a matrix 
entry which could be zero or non-zero is denoted by a 
star (+).] 


‘Type of solution 
ve ele 
Ov «le Unique solution 
00 vi« 
No zero row of new (case (i)) 
coefficient matrix 
[: Ys ‘| No solution 
0 0 olv 


A zero row of new coefficient | (equations inconsistent; 


matrix, but not of new case (ii)) 

augmented matrix 

(i.e. Equation Esp has form 

0 = non-zero number) 
Ya s/s | ae 
ueiice wells An infinite number of 
0 0 olo solutions 

‘A zero row of the new (case (iii)) 


augmented matrix 
(ie. Equation Ey, has form 
0=0) 


In case (iii) we can obtain the solution set by taking 
3 = k, where k is an arbitrary number, then 
back-substituting in Equations E2, and E; as usual. 
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2. The rows R1,R2,...,R» of a matrix are said to be 
linearly dependent if some linear combination of 
these rows is zero, i.e. if constants ki,k2....,kn (not all 
zero) exist such that 


ky Ry + koRe +--+: + kpRn = 0, 
where O here means a row of zeros. 


If the rows of the matrix A of coefficients in the given set 
of equations are linearly independent, then the solution is 
unique (case (i)). If the rows of A are linearly dependent 
but those of the augmented matrix A[b are linearly 
independent, then there is no solution (case (ii). If the 
rows of Alb (and hence also those of A) are linearly 
dependent, then there is an infinite number of solutions 
(case ({ii)). 


Section 3: Some possible difficulties 


1. (i) The kth pivot (or pivot element) is the one by 
which we divide to obtain the multipliers in the Ath 
step of the elimination. 


(ii) ‘The pivot row is the row in which the current 
pivot lies. 


(iii) Essential row interchange is the process of 
interchanging the rows of a matrix when the 
obvious pivot is zero, This can be done by 
interchanging the current pivot row with whichever 
of the rows beneath it will provide the pivot with 
the largest modulus. 


(iv) Partial pivoting is the process described in the 
preceding sentence with the difference that it is 
done even if the obvious pivot is not zero, provided 
that it will replace the obvious pivot by a pivot 
with a larger modulus. It ensures that no multiplier 
has modulus exceeding 1. 


2. The section shows some of the difficulties presented by 
simultaneous equations and their solution. 


(i) Il-conditioning. A set of equations is said to be 
ill-conditioned if a ‘small’ change in the data (the 
coefficients, or the right-hand sides of the 
equations) produces ‘large’ changes in the solution. 
Il-conditioning is a property of the equations 
themselves, not of the method used to solve them. 
One way to detect ill-conditioning is to vary the 
coefficients slightly and see if there is much change 
in the solution. A possible cure for ill-conditioning 
is to replace one of the equations by an equation 
containing some new independent information, and 
hence obtain a set of equations which are not 
ill-conditioned. 


(ii) Induced instability. This occurs when we start 
with a set of equations which are not ill-conditioned 
and yet the method used to solve the equations 
induces ill-conditioning in the problem during the 
process of solution. It is usually heralded by the 
appearance of a very large multiplier. Partial 
pivoting is used as a way of trying to avoid this 
(and to avoid the consequent build-up of error). 


Section 4: The computer package 


Section 5: Special cases 


If the right-hand sides of the equations are all zero, the 
equations are said to be homogeneous. If the number of 
unknowns is the same as the number of equations, and 
the equations are linearly independent, then the only 
solution is the trivial solution, in which all unknowns 
equal zero. If, however, the equations are linearly 
dependent, then there is an infinite number of solutions, 


A set of m equations in n unknowns has an infinite 
number of solutions if m <n (unless the equations are 
inconsistent) and no solutions at all if m > n (unless at 
least m —n of the equations are linearly dependent on 
the others). 


Unit 12 Heat transfer 
Section 1: The nature of heat transfer 


1. Heat transfer means energy transfer by heating 
and/or cooling. This kind of energy transfer requires the 
existence of a temperature difference, and the energy is 
transferred from regions or bodies of higher temperature 
to those of lower temperature. There are three modes of 
heat transfer: conduction within a (normally stationary) 
material, radiation (or electromagnetic radiation) which 
does not require any material, and convection by 
motion of a fluid. In free convection the fluid motions, 
or convection currents, are driven by density changes 
caused by the temperature difference itself; in forced 
convection, they are driven by some other agent. 


2. Part of the energy of a body depends on its 
temperature. This part is known as its internal energy, 
and the change of internal energy due to a change from a 
temperature 6; to a temperature 2 is me(@2 — 61), where 
m = mass of body, 
= specific energy capacity 
or specific heat of the material. 
Power is the rate at which energy (of some particular 
kind) is produced. 


Section 2: Steady-state conduction in one 
dimension 


1. Fourier’s Law for steady-state conduction in one 
dimension is 
do 
=-KAS, 
q=-KAT 
where 


q = heat transfer rate in the positive x direction, 

« = thermal conductivity of the material, 

A = area at right angles to the direction of heat transfer, 

6 = temperature, 

x = distance measured in direction of temperature 
variation. 


4 is known as the temperature gradient. 


Steady-state means that the temperature at any given 
point does not change with time; one-dimensional 
means that it depends on only one of the three possible 
space coordinates (in a suitable coordinate system). 


2. For steady-state conduction through a uniform slab of 
constant cross-sectional area A, Fourier's Law reduces to 


A — 

— 
where 6, — 02 (for 0, > 62) is the temperature drop 
across thickness b of the material. 


q=KA 


For steady-state conduction through a pipe of circular 
cross-section with internal radius r;, external radius r2 
and length /, Fourier’s Law leads to 

_ wns — 62) 

~ Tog.(r2/r1) * 
where ; — 62 (for 6; > 02) is the temperature drop 
across the pipe wall. 


Section 3: Walls and radiation 


1, Heat transfer through walls and roofs is often 

represented by steady-state one-dimensional models. 

The U value of a wall is defined as follows: 

rate of heat transfer per unit area 
overall temperature difference 

New techniques and components are being developed to 

satisfy the demand for very low U values in buildings. 


U value = 


2. The thermal energy emitted by radiation from a body 
depends on the fourth power of the absolute 
temperature of the body, ie. 
thermal energy emitted by radiation = 6, 
where 
jis a constant, 
0, = 273.2 + @ = absolute temperature (in kelvins), 
@ = temperature in °C. 


Section 4: Convection and insulation 


1, Convective heat transfer in a fluid can be modelled 
using the following assumptions: 


(i) the temperature in the bulk fluid is uniform, say 6; 


(ii) if the fluid meets a plane solid surface with area A 
and temperature @,, then the rate of heat transfer 
from surface to fluid is hA(@, — @;), where h is the 
convective heat transfer coefficient. 


2. The equation for heat transfer through a wall or 
window with area A is 

q = AU (An — Pour), 
where On and Aour (for Ain > @our) are the bulk air 
temperatures inside and outside the building. Overall 
U values which take both conduction and convection into 
account can be calculated by adding formulas for 
adjacent temperature differences. 


In particular, for a solid wall 
eet eran day ba 
alee 
where b and x are the thickness and thermal conductivity 
of the wall and hjy and hou are the convective heat 
transfer coefficients at the inside and outside surfaces of 
the wall respectively. 
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3. Convective heat transfer from a surface such as the 
wall of a pipe can be increased by means of fins. The 
heat transfer in the fin can be modelled by assuming that 
the temperature in the fin depends only on its distance 
from the pipe wall, dividing the fin into elements, and 
then considering convection from each element of the fin 
to the surrounding air and conduction from each element 
to its neighbours. This leads to the differential equation 

@@ 2h 

Pr ti Pa — 8a), 
where h is the convective heat transfer coefficient, « is 
the thermal conductivity of the material of the fin, b is 
the fin's thickness and 6, is the temperature of the air. 


4. The heat transfer properties of a layer of fluid (such as 
the air space in double glazing) are complicated, and are 
best modelled using a combined heat transfer 
coefficient h., which is defined analogously to a U value 
but depends on the temperature difference across the 
layer as well as on its thickness. 


Unit 14 Vector algebra 
Section 1: Vectors: scaling and addition 


1. A (non-zero) vector is a mathematical object 

consisting of two parts: 

(i) @ positive real number called the modulus or 
magnitude of the vector, and 


(ii) a direction in space (or in the plane). 


The modulus of a vector v is denoted by |v| or sometimes 
by uv. We also define the zero vector to have modulus 
zero; no direction is defined for it. The zero vector is 
denoted by 0. In the context of vectors, we use the word 
scalar to denote a real number. 


2. Physical quantities that can be modelled by vectors or 
by scalars are called vector quantities or scalar 
quantities respectively. Examples of vector quantities 
include displacements (i.e. changes in position), velocities 
and forces. The displacement from a point P to a point 


Q is represented by the displacement vector PQ. 


3. Equal vectors have the same magnitude and direction. 
Scaling of a vector is multiplying it by a positive 
number; this leaves the direction unaltered and multiplies 
the magnitude by that number. Multiplying any vector 
by the number zero gives the zero vector. Multiplying a 
vector v by a negative number —m gives a vector —mv 
with magnitude m|v| and direction opposite to that of v. 
This new direction has the same orientation (i.e. in two 
dimensions, the same slope) as that of v but its sense is 
opposite to that of v. In particular the vector —v has the 
same magnitude as v and the opposite direction. 


4. A unit vector is one whose magnitude is 1. The 
vector at is a unit vector with the same direction as a. 
Unit vectors along the Cartesian axes are called 
Cartesian unit vectors, and are denoted by i, j, k. 
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5. Vector addition is defined by a parallelogram rule 


or alternatively by a triangle rule 


Vector subtraction is defined by 
a—b=a+(-b). 


6. Vector addition is commutative and associative (see 
Subsection 3.6 of this Handbook), and distributive with 
respect to multiplication by a sealar: 

(m, + m2 )a = ma + maa, 

m(a +b) = ma+ mb. 
Multiplication by a scalar has the ‘associative’ property 
that 

m(maa) = (mymo)a. 
‘The zero vector has the properties 0a = 0 and 0+a=a. 


Section 2: Components of a vector 


1. Given a Cartesian coordinate system, any 
two-dimensional vector a can be written uniquely in the 
component form 

a=ait+aaj, 
where i and j are unit vectors along the positive z- and 
y-axes respectively. The vectors ai and aaj are called the 
(Cartesian) component vectors of a, and the numbers 
a; and a» are called the (Cartesian) components of a 
in the directions of i and j, The magnitude of a is 
s/a? + a3. To add two vectors, we add the respective 
components; to multiply a vector by a scalar, we multiply 
the respective components by that scalar. 


2. For three-dimensional vectors the representation 
a=aji+a2j+ask 
has analogous properties, where k is a unit vector along 
the positive z-axis, chosen at right angles to the x- and 
y-axes. A right-handed coordinate system is one 
satisfying the right-hand rule: if the thumb of your 
right hand points along i and the forefinger along j, then 
the middle finger can comfortably point along k. A 
right-handed coordinate system also satisfies the 
right-hand screw rule that a screw rotated from i to j 
will advance in the direction of k. 


3. Ifa point P has coordinates x, y, = relative to a given 
Cartesian coordinate system with origin O, the 
displacement vector OP is called the position vector of 
P, and is given by 

r=zi+yj+2k. 


Section 3: Products of vectors 


1. The dot product or scalar product of two vectors 
a and b is the scalar defined by 

a+ b= |al [b|cos8, 
where (for 0 < @ < 7m) is the angle between the 
directions of a and b. 
In particular we have 

a-a=|al’. 
The angle @ between two non-zero vectors a and b is 
given by 


2. The projection of a vector b onto a given direction is 
defined to be bj cos, where 0 is the angle between the 
direction of b and the given direction; it is equal to 
a-b/jal, where a is any vector whose direction is the 
given direction (it is often convenient to make a a unit 
vector). Hence 


a+b = [al times the projection of b onto the direction of a 
= |b| times the projection of a onto the direction of b, 


3. The dot product is commutative, 
a-b=b-a, 

and distributive, 
a:(b+c)=a-b+a-c, 
{a+b)-c=a-c+b<c, 


but it is non-associative and there is no inverse. It 
combines with sealing according to the rule 


(ma) +b = m(a+b) = a+ (mb). 
If a and b are non-zero vectors, then a+b = 0 if and only 
if a and b are perpendicular vectors. 


4. Dot products of Cartesian unit vectors are 
si=jej=k-k=1, 
+j=j-k=k-i=0. 
If two vectors a and b have Cartesian components 
@},a3,a3 and b;,b2,by respectively, then 

a-b= ah; + azby + abs. 


The formula |a/? = a+ a= a? + a3 + a3 is a particular 
case of this result. 


5. The cross product or vector product of two 
vectors a and b is the vector defined by 


ax b= [al |bisin@c, 


where @ (for 0 < @ < 7) is the angle between the 
directions of a and b, and c is a unit vector at right 
angles to both a and b. The sense of c is given by the 
right-hand screw rule, that a screw rotated from a to 
b advances along the direction of c. 


6. The cross product is distributive, 
ax(b+c)=axb+axc, 
(a+b)xc=axc+bxe, 

but is non-commutative, obeying instead the law 
axb=-(bxa). 


It is non-associative, and there is no inverse. It combines 
with scaling according to the rule 

(ma) x b = m(a x b) = a x (mb). 
If a and b are non-zero vectors, then a x b = 0 if and 
only if a and b have the same orientation (i.e. they are 


parallel or anti-parallel). In particular, we have 
axa=0. 


T. Cross products of the Cartesian unit vectors are 
=kxk=0, 
jxk=i, kxi=j. 
If two vectors a and b have Cartesian components 
1,42, a3 and by,b2,bs respectively, then 
a x b = (a2b3 — asba)i + (asbi — a1bs)j + (aib2 — aabi )k. 


8 The triple scalar product of three vectors a, b, ¢ is 
a-(b xc), 


9. Given a plane, any vector n perpendicular to it is 
called a normal vector, and the equation of the plane 
can be written 

(w — a)ny + (y — b)na + (2 —€)ng = 0, 
where (a,6,¢) is a point in the plane and 
n=mi+naj+nsk. If a and b are vectors lying in this 
plane, with different orientations, then a x b is a normal 
vector. 


Section 4: Differentiating vectors 


1, A vector function is a function whose domain is the 
set of all real numbers and whose codomain is the set of 
all vectors: to each real number ¢ the vector function f 
assigns a vector f(t). Alternatively, the formula for f{t) 
may itself be described as a (vector) function of t. The 
Cartesian components of a vector function are real 
functions of the usual kind: 


f(t) = Althi+ fA(Oj + flOk. 
‘The magnitude and direction of f(t) are also functions 
of t. A vector whose Cartesian components are constants 
(and hence whose magnitude and direction are constants) 
is called a constant vector. 


2. As t varies, the point whose position vector (relative 
to a given coordinate system) is f{t) traces out a curve in 
three dimensions, which helps in visualizing the vector 
function f. 


3. The derivative of a vector function f is 


df(t . t+h)—fit 

a= ag ( E22 -©). 

If f(t) has the Cartesian component form 
F(t) = fi(ti+ fol) + fo(tk. 

then 

The rules for differentiating sums and products are 


if eartiviieset a 8" 
dt dt 
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diay = vou, 
feuxv= (xv) + (ax), 


4. The tangent at any given point on the curve 
representing a vector function f is in the direction of the 


df(t) 


vector —— at that point. If r(t) is the position vector of 
dt 
dr(t) 
.- Ifa 
vector function g(t) has constant magnitude, then 


dg _ 
eee 


a particle at time t, then its velocity vector is 


Unit 15 Newtonian mechanics in 
three dimensions 


Section 1: The ingredients of Newton’s second 
law 


1. In three-dimensional Newtonian mechanics the 
position of a particle is represented by the position 
vector 


r=zi+yj+zk, 
where x, y and = are the coordinates of the particle, and 
i, j and k are unit vectors along the axes of a static 
right-handed Cartesian coordinate system that is 
calibrated in metres, 


2. The velocity vector v is found by differentiating the 
position vector with respect to t, the time measured in 
seconds: 


v=t=ai+ pj+ik. 
‘The magnitude of the velocity vector is the speed of the 
particle. The direction of the velocity vector is the 


direction of motion of the particle. 


3. The acceleration vector a is found by 
differentiating the velocity vector with respect to t: 

= #i+ jj+ ik. 

This describes the rate of change of velocity both in 
magnitude and direction. The acceleration will be 
non-zero if either the speed or the direction of motion 
changes. 


4, Bach individual force acting on a particle is 
represented by a force vector. If the direction of the 
force makes angles a, 3 and with the positive «-, y- and 
z-axes, and if the magnitude of the force is |F|, then the 
force vector is 

F = |F|(cos ai + cos 8j +.cosk). 


The cosines that appear in this equation can sometimes 
be simplified by using the trigonometric identities 


cos (Z—4) =sing, cos (+0) = sing, 
cos(x — ¢) = —cos¢. 
5. The gravitational force vector for a particle of 


mass m is mge, where e is a unit vector pointing 
vertically downwards and g ~ 9.81ms~*. 
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6. Ifa particle experiences simultaneously a number of 
individual forces, represented by the vectors F,F2,..., it 
behaves just as if it were acted upon by the single force 

F=Fi+Fo+---. 
This formula is known as the law of addition of forces, 
and the vector sum F is known as the total force 
vector, The components of the total force vector are 
denoted by F;, Fy and F:. 


Section 2: Newton’s second law 


1. Newton’s second law, in its vector formulation, 
states that, at each instant of time, the motion of a 
particle is governed by the equation 


= re dr 
mi=F, ie. maa =F 
or, in terms of components, 
mé=Fe, mi = Fy, 


where 
m is the mass of the particle, 
¥ = #i + jj + Zk is the acceleration vector, 
F = F,i+ F,j + F.k is the total force vector. 


2. A particle that moves in a circle with constant speed 
is said to perform uniform circular motion, and its 
acceleration is directed towards the centre of the circle. 
For example, a particle of mass m for which 

r(t) = U(coswti+ sinwt j), 
where w is a constant, performs uniform circular motion 
at speed lw round a circle of radius /. Its acceleration is 
—w'r(t) and the force on it is —mw*r(t), which has 
magnitude mw?! and is directed from the particle towards 
the centre of the circle. 


Section 3: Algebraic equations and sloping 
tables 


1. A static particle, i.e, one which is at rest over some 
interval of time, experiences zero total force. 


2. A particle in contact with a solid surface experiences a 
normal reaction force, Fy, and a force due to friction, Fy. 


3. The normal reaction rule: Fy is perpendicular to 
the surface and |Fy| is just sufficient to prevent the 
particle from sinking into the surface. For the particle to 
remain in contact with the surface, |Fy| > 0. 


4. The force of friction, or frictional reaction, Fy, is 
parallel to the surface, in the direction that opposes 
motion. For a static particle, |F,| is just sufficient to 
prevent sliding over the surface, but the particle can 
remain static only so long as |Fy| < |F |, where p is the 
coefficient of static friction. For a moving particle, 
|Fy| ='|Fy|, where y’ is the coefficient of kinetic 
friction. 


Section 4: Uncoupled differential equations 
and shot-putters 


1. We can model a projectile, such as the shot, as a 
particle that experiences only the downward force of 
gravity. The path or trajectory of a moving point is the 


set of points it passes through. The path of the shot is a 
parabola, that is, the graph of a quadratic function 
(with the y-axis vertical). 


2. If the shot-putter launches from zero height, the 
maximum range is 
Vinax 
nae 
where vmax is the greatest speed with which he can 
release the shot. ‘The maximum range is achieved for an 
angle of release of }7 radians. 


3. If the shot-putter launches from a height h, the 
maximum range is 


Ly1+2h]L, 


where L = v3iax/g, and this is achieved for an angle of 
release of 


ee 


Section 5: Coupled differential equations and 
pendulums 


1. This section discusses a particle of mass m suspended 
from a fixed point by a light taut inextensible string, 
where ‘light’ means ‘of negligible mass’, ‘inextensible’ 
means that l, the length of the string, is constant, and 
‘taut’ means that the distance between the two ends of 
the string is /. 


2. If the particle moves in a horizontal circle, then its 
speed is constant, the tension in the string is constant, 
and the equations of motion for x, y and z are uncoupled 
(i.e. the equation for # does not involve y or 2, and so on). 


3. If the particle moves in a vertical circle, then its speed 
is not constant, the tension in the string is not constant 
and the equations of motion are coupled. In this case the 
angle @ between the string and the downward vertical 
satisfies the equation of motion 


@ 


$ sind. 

This is consistent with the energy equation 
}m(16)* + mgl(1 — cos 0) = E = constant. 

The tension in the string is 


T= ml + mgcosd 

= mg(3cos 0 — 2) + 28 
so long as the right-hand side of this equation is positive. 
As soon as T = 0 the string goes slack, and from then 
onwards a different description of the motion must be 
used. 


4. The sum of the kinetic and gravitational potential 
energies of a particle is conserved provided the particle 
experiences only the force of gravity and forces (like the 
tension in a string, the normal reaction force or the 
magnetic force on a charged particle) that act at right 
angles to its velocity. 


Unit 17 The dynamics of 
many-particle systems 


Section 1: Two-particle string—pulley systems 
Section 2: Systems of particles 


1, Newton's third law states that the inter-particle 
forces acting mutually between a pair of particles are 
equal in magnitude and opposite in direction. If the 
particles are labelled i and j, then this law may be 
expressed in terms of vectors by the equation 

Fy = -Fyi, 
where Fy is the force on Particle i due to the effect of 
Particle j, and vice versa. 


2. For a system of n particles with masses 
m1, M2,...,My and respective position vectors 
T1,T2,-+-)Pn, the centre of mass has position vector 


1 x 
R= women + mare +-+++ Maran) = M dre 
where 
M =m +ma+--+my = Som 
‘=I 
is the total mass of the system. If the particles are acted 
upon by external forces F{*', FS**,..., F&** respectively, 


and if the inter-particle forces satisfy Newton's third law, 
then the equation of motion of the centre of mass is 


eer et cere = See 


is the total external force acting on the system. 


3. In the absence of external forces, the equation of 
motion of the centre of mass reduces to R = 0, indicating 
a uniform motion 

R(t) = Vot + Ro, 
where the constant vectors Vo and Ro are respectiv. 
the velocity of the centre of mass and its initial position 
(at time t = 0). 


4. The centre of mass of a symmetric homogeneous 
object is located at the geometric centre of the object. 
The centre of mass of a collection of objects can be found 
by treating each object as if it were a single particle 
situated at its own centre of mass. 


Section 3: Linear momentum and energy 


1. The linear momentum p of a particle with mass m 
and position vector r is defined by 

p= mi. 
In terms of linear momentum, Newton’s second law 
can be expressed as 

p=F, 
where F is the total force acting upon the particle. In the 


35 


absence of any force on the particle, its linear momentum 
is conserved (that is, p is constant), If the mass m of the 
particle is constant, then Newton's second law reduces to 


mt =F. 


2. For a many-particle system, the total linear 
momentum P is the sum of the momenta of the 
individual particles. It is related to the total mass M of 
the system and the centre of mass R by the equation 
P=MR. 
Newton's second law leads to the result 
P=F™ 
where F® js the total external force acting on the 
system. In the absence of external forces, the total linear 
momentum is conserved (that is, P is constant). If the 
total mass M of the system is constant, then the last 
equation reduces to 


MR=F™. 


3. For a system of particles involved in a collision in the 
absence of external forces, momentum conservation 
implies that Ppefore = Pater. In the two-particle case, 
this is equivalent to 

my Uy + mzuz = m1 Vi + M2v2, 
where the particles have masses m1, mo, initial velocities 
uy, up and final velocities v1, v2, respectively. 


4. Ifa two-particle collision is elastic, then the kinetic 
energy equation 

4m, us|? + dma|ua|? = dmalvil? + malval? 
is also satisfied. Otherwise the collision is inelastic. 


5. Newton's law of restitution for a two-particle 
collision is 

(v2 — v1) -1=—e(u2 — m1) +1, 
where 1 is a unit vector in the direction of individual 
particle momentum change, and the constant e is the 
coefficient of restitution for the objects represented by 
the particles. For an elastic collision, the coefficient of 
restitution is e = 1. If the collision is one-dimensional, 
then this law reduces to 


v2 — v1 = —e(ua — ur). 


Section 4: Newton’s third law 
Section 5: Rocket motion 


1. The velocity v of a rocket satisfies the rocket 
equation 

mv —mu=F™, 
where the rocket (with unspent fuel) has mass m, is acted 


upon by an external force F*** and emits exhaust gases 
with relative velocity u. 


2. In the absence of any external force, the rocket 
equation reduces to 


ae 
v=—u. 
m 
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Unit 18 Polynomial approximations 


Section 1: Basic theory 


1. A polynomial of degree n is a function that can be 
expressed in the form 

dnt" + ane"! +--+ + aie +o, 
where a, #0. The constant numbers ay, ay—1,-.- 
are called the coefficients. 


31,09 


2. To evaluate the above polynomial for a given value of 
x by nested multiplication, use the recurrence relation 


Urs) = Up +A@n—r-1 fore =0,1,2,....2—1, 
with uo =an. Then tn is the value of the polynomial. 


3. For a given function f(x) and a given number a in its 
domain, the nth Taylor polynomial for f about a is 
1 ” 
Sa) + (@—a)f'(a) + (2 - @)*s"(a) 


+ A (e=a)" f(a) = ys Ae- ay fa). 


4. Taylor’s Theorem. Let the function f and its first 
n+ 1 derivatives exist and be continuous at all points 
between and including two given numbers a and x. Then 
we have 

L(x) = p(x) + r(x), 
where 


we) = > leah 7a) 


is the nth Taylor polynomial for f about a, and r(x), the 
remainder, satisfies 


(e- a)"tt ft(c,) 


for some number cx between a and x. Taylor's Theorem 
with n = 0 is also called the Mean Value Theorem. 


5. The error function for the approximation of a given 
function f(2) by a polynomial p(x) is 

e(x) = p(x) — f(z). 
In the case of Taylor approximations it is related to the 
remainder by 


e(x) = —r(z). 


6. Suppose f is a given function and a,b, M are numbers 
such that 

If" (c)| <M 
for all c satisfying a <c <b. Then, for any a satisfying 
a<a<b, the error in the nth Taylor polynomial 
approximation to f about a satisfies the error bound 
formula 


le(x)| < wim -al"'M fora<ez<b. 


Section 2: Two applications of Taylor 
polynomials 


1. A root of the equation f(z) = 0 means a solution of 


that equation. 


2. To calculate a root of the equation f(x) = 0 by the 
Newton-Raphson method, use the recurrence relation 
S(zr) 

Fee)" 

Tf the initial guess zy is not too far out, the sequence 
©,21,£2,... will converge to the root (i.e. for sufficiently 
large r, the value of z, is arbitrarily close to the root). 
The error in x,+1 is proportional to the square of the 
error in z,; so if the sequence converges at all it converges 
rapidly. 


Try = Pr — 


3. The Taylor series method of order n for the 
solution of the first-order differential equation 
y =m(z,y); 
where the value of y(xo) is given, involves the following 
steps. 
(i) Differentiate the expression m(z,y) with respect to 
x (remembering that y is a function of x) n — 1 
times to obtain expressions for 


di 
law) fm y), co ey). 


os 
(ii) Choose a step size h. 
(iii) Calculate Y;, Ya,... 


Your = Yet hYe + MY! + 


from the recurrence relation 

va Spry, 
nt 

where 


tr =zo+Th, 
ee peas ds 


Y= om 2..¥e)s 


yin) = = (2, Yr); 


with the initial condition 
Yo = y(z0). 
(iv) ¥; is an approximation to y(z,). 


Section 3: The catenary 


1. Hyperbolic cosines and sines are functions defined 
by 

cosh x = $(e* +e™*), 

sinh z = 3(e* — e~*). 
‘They have the properties 

cosh? x —sinh* x = 1, 


A cos 2 = sinh, 


& sinh 2 = cosh 2, 


2 
coshe = 145+ 5 


P z 
sinhe=2+3, + 


‘These functions are related to the hyperbola x? — y’ 
in the same way that the ordinary cosine and sine 
functions are related to the circle. 


2. The length L of a curve with equation y = y(z), 
between two points on it with x-coordinates a and , is 


ux fas () ee. 


3. The curve with equation 
1 
= =coshar, 
y= {cos 


where a is any positive constant, is called a catenary. It 
models the shape of a rope suspended between two points. 


Section 4: Interpolation and integration 


1. If the values of a function f(z) are known at points 
@o,%1,.+.,in, arranged in increasing order but not 
necessarily equally spaced, then the use of some other 
function g(a), which takes the same values as f(x) for 
© = Z0,21,.--,2n, to approximate f(x) for values of x 
between xp and zn, is called interpolation; the 
corresponding procedure for values of x less than xo or 
greater than x, is called extrapolation. 


2. There is a unique interpolating polynomial of 
dogree <n which takes the same values as f(x) for 
© = 10,01,...,2n- Its n+ 1 coefficients, say 

49, @1,...,n, can be found by solving the n +1 
simultaneous linear equations 


ao + ara +++: +anz" = fli) (i =0,1,...,n). 


8. Provided xo < x < ap, the error function e(x) for the 
approximation of f(a) by the interpolating polynomial 
defined above satisfies 


e(x) = — (w — xo)(x — 21)...(2— an) f"* (ce), 


a 
(n+1)! 
where cs is some number satisfying xo < es San. If M is 
a number such that 


[f*P(c)| <M 


for all c between xo and x», then e(x) satisfies the error 
bound formula 


lea) < anil = 20)(e — 21)... (@ —24)|M. 


4. Buler’s integration formula, based on 
approximating the area under a segment of a curve by a 
rectangle, is 


0 
if f(a) dx = h(f (wo) + f(e1) +++ + L(en-1)), 


where 
h=(b—a)/n 

and io,21,...,2tn—1 are equally-spaced points defined by 
2 =atrh (r=0,1,....n—1). 


5. The trapezoidal method for integrals, based on 


approximating the area under a segment of a curve by a 
trapezium, is 


» 
r) F(a) dx ~ h(5 f(20) + f(e1) +--+ + f(@n—1) + 3S (2n)), 


where 
h=(b~a)/n, 
t=atrh (r=0, 


yn). 
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6. Simpson’s method for integrals, based on 
approximating the area under a segment of a curve by a 
parabola, is 


: 
[ $(a)dz ~ 3h(f(20) + 4f (zs) + 2f a2) + 4f(2s) 


+ 2f (za) +++ +4F(tn-1) + f(2n)), 
where 
h=(b—a)/n, 
z=a+rh (r=0,1,...,n), 
and n must be even. 


7. The three integration formulas are examples of 
Newton-Cotes formulas, obtained by integrating an 
interpolating polynomial. They are said to be in 
composite form if n > 1 for Euler's and the trapezoidal 
formulas, and if n > 2 for Simpson’s. For a given step 
size h, Euler's formula is normally the least accurate of 
the three and Simpson's the most accurate. 


Unit 19 Numerical methods for 
differential equations 


In all these methods the differential equation is taken to 
be 

y' =m(z,y), 
with y(xo) given, and we seek approximations to y(a) for 
@ = 21,22,..., where 

t =zo+rh 
and h is the step length or step size. The approximation 
to yr, which means y(z,), is denoted by Y,. The 
approximation to y/, which means y/(,), is denoted by 
Y/ and is given by the formula 


Y¥; = m(z,,¥-). 
Section 1: Numerical methods 
1. By an integration method for the numerical 
solution of differential equations we mean a method 


based on a numerical integration formula. 


2, Buler’s method uses the recurrence relation 
Yeu = ¥ + hz. 


3. The trapezoidal method uses the formula 

Your = ¥, + $h(¥y + ¥e41)- 
In the case of a linear differential equation, for which 

m(z,y) = U(z)y + k(x), 
the formula of the trapezoidal method can be solved for 
Y-41 to give 

(1+ EAl,)¥, + Sh (hy + hepa) 
Yoo A 
Nya 

where I, =I(z-), kr = k(z-), and so on. 


4, An explicit method is one where the recurrence 
relation for ¥-41 does not involve ¥/,1, ¥41, or higher 
derivatives at x,41. An implicit method is one where 
the recurrence relation does involve Y;/,, or higher 
derivatives at 2-41. 
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5. The local truncation error is the error introduced 
in a single step, estimated as the value of ¥-+1 — yr+i 
when Y, = yr. The principal term in the local 
truncation error is the first non-vanishing term of the 
Taylor series about x, in powers of h for the local 
truncation error; this series is calculated by subtracting 
the Taylor series 


Vet = Ye thy, + Shy! +. 
from the corresponding series for ¥;-+1 obtained from the 
recurrence relation. 


6. The principal term in the local truncation error for 
Euler's method is —}h?y’, for the Taylor series method 


of order n (with n > 2) it is —A"*1y\"*" /(n + 1)!, and 
for the trapezoidal method it is 5h°yy". 


Section 2: The predictor-corrector method 


1. The Euler-trapezoidal predictor-corrector 
method calculates Y,+1 in two stages: 


m1 = Ye + hYy where Y/ = m(z,,Y-), 


You = Yet SA(¥/ + ¥eth1) where ¥ef1 = m(ze41,¥eis) 


The principal term in the local truncation error for this 
method involves h*. 


2. In general a predictor-corrector method is one 
where we calculate ¥,41, the approximation to yr+1, in 
two stages: the first stage uses an explicit method to give 
a first approximation Y,",,, which is then used on the 
right-hand side in the recurrence relation of a more 
accurate implicit method to give an improved 
approximation Y;1. 


Section 3: The analysis of numerical methods 


1. A one-step numerical method for the solution of 
the differential equation y’ = m(z,y) is a method which 
approximates the differential equation by a first-order 
recurrence relation of the form 
Yeti = ¥r +h O(ar, Yrs Yeo, h). 
It is consistent with the differential equation if 
(Er Yrs Yrs 0) = m(zr,Hr)- 
An alternative statement of the definition of consistency 
is 
. Yeti = Yr 
lim h 


jim = m(zr.¥r)- 


2, For a fixed value of x, say x*, the.global error of a 
numerical method for calculating y(x*) using N steps of 
length h = (x* — x)/N is defined as Yx — y(x"). 

The method is said to be convergent on some interval 
{zo, 6] if for all x* in this interval (i.e. all x” satisfying 
xo <2" <b) we have 


jim Yn = y(x"). 


3. A one-step method is convergent if and only if it is 
consistent. Further, if the principal term in the local 
truncation error is h?*", for some integer p, then the 
global error at 2” is given approximately by 

Yn —y(a") = Ch?, 
where, for small h, C does not depend on h. 


In particular, the global error for Euler's method is 
approximately proportional to h, whereas the global 
errors for the trapezoidal method, the Buler-trapezoidal 
method and the Taylor series method of order 2 are all 
approximately proportional to h? for small h. 


4. A one-step method, applied to a given differential 
equation problem, is said to be absolutely unstable if 
the resulting recurrence relation problem is absolutely 
ill-conditioned, and to be absolutely stable if the 
recurrence relation problem is absolutely 
well-conditioned, 


5. We can test the stability of a given numerical method 
by applying it to the test problem y/ = ay with 

y(zo) = 1, where a is a constant whose value is to be 
chosen. The interval of absolute stability of the 
method is the set of values of ha for which the method is 
absolutely stable. 


6. The interval of absolute stability for Euler's method, 
for the Euler-trapezoidal predictor-corrector method and 
for the Taylor series method of order 2 is (—2,0), i.e. 
these methods are stable if -2 < ha < 0. The interval of 
absolute stability for the trapezoidal method is (—90,0), 
i.e, the method is stable for all negative values of ha. 


7. Applied to the equation y' = m(z,y), a given 
numerical method is absolutely stable if h™ (a, Ye) lies 


in the interval of absolute stability for all values of r used 
in the calculation. 


Section 4: Simpson’s method 


1. Simpson’s method uses the recurrence relation 
Your = Yeo + $h(¥yoy + 4Yy + Yrs): 

Since the recurrence relation is of second order, this may 

be called a two-step method. 


2. Simpson’s method is implicit, but if it is applied to a 
linear differential equation 
y =U(x)y + k(x), 
the recurrence relation can be solved for Y¥-+1 to give 
~All Ye + (3 + lt) ¥r—1 + (hr + he + Bega) 
a 3— Alaa . 
where L, = l(a), ke = k(2r), and so on. 


3. The principal term in the local truncation error for 
Simpson's method contains the factor h*, and the method 
is therefore capable of very high accuracy. 


4. A spurious solution of the recurrence relation for 
Simpson's method is a term in the general solution which 
does not correspond to any solution of the differential 
equation. For the differential equation 


y =Ux)y + kz), 
the spurious solution is harmless provided I(x) is positive, 
but if I(x) is negative the spurious solution makes the 
method unstable. 


Unit 20 Matrix algebra and 
determinants 


Section 1: Some simple matrix operations 


1. Anm xn matrix is a rectangular array of numbers 
called elements, in m rows and n columns. The two 
numbers m and n together specify the size of the matrix. 
The element in the ith row and jth column of the matrix 
A can be denoted by ajj. A column matrix or column 
vector is a matrix with just one column; a row matrix 
or row vector has just one row. 


2. Two matrices A and B are equal if they are the same 
size and corresponding elements in each matrix are equal. 


3. If A is a matrix and k is a real number, then kA is 
the matrix of the same size whose elements are k times 
the corresponding elements of A. 


4, If two matrices A and B are the same size, then we 
can add them, the sum A + B being the matrix formed 
by adding the corresponding elements of A and B. 


5. A matrix whose elements are all zero is called a zero 

matrix and is written 0. For any matrix A we have 
A+0=0+A=A 

and 0A=0. 


6. Given any matrix A, the matrix (—1)A, whose 
elements are the negatives of those of A, is written —A. 
‘This matrix has the property 
A+(-A) =0, 
and it can be used to define matrix subtraction as 
A-B=A+(-B). 


7. The algebraic operation of adding matrices of a given 
size is closed, commutative and associative, and is 
distributive with respect to multiplication by a scalar (see 
Section 3.5 of this Handbook). 
Multiplication by a scalar obeys the further rules 
(kike)A = ki(k2A), 
(ki + k2)A = ky A + ke A. 


Section 2: Matrix multiplication 


1. If A is an m x n matrix and B is an n x p matrix, 
then we can multiply them, with A written on the left. 
The product C = AB is an m x p matrix whose element 
in the ith row and jth column is 


2. Matrix multiplication is associative, and is distributive 
with respect to matrix addition, i.e. 


A(BC) = (AB)C, 
A(B+C)=AB+AC 
and (A+B)C=AC+BC. 
However, it is not (in general) commutative: 
BA # AB. 
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3. The product of a matrix and a column vector, with 
the matrix on the left, is another column vector. A set of 
simultaneous linear algebraic equations can be written in 
the form 

Ax =b. 
where A is the (given) matrix of coefficients, x is a 
column vector whose elements are the unknowns, and b is 
a given column vector. Moreover, if A is a given matrix 
and x is a variable column vector of suitable size, then we 
can say that A transforms x to Ax (i.e. A defines a 
function whose domain is a set of column vectors; this 
function is called a linear transformation). 


4. The transpose of a matrix A, written A’, is 
obtained by interchanging the rows and columns of A, 
ive. the element in the ith row and jth column of AT is 
a,;. The main properties of transposed matrices are 
(AT)T=A, 
(AiA2,..An—A 


Section 3: Change of axes 


Suppose Ox, Oy and O'x', Oy’ are two Cartesian 
coordinate systems in the plane, with origins O and O’, 
that the origin O' of the second system has coordinates 
(a,b) with respect to the first system and that the 
direction of the O'z' axis is obtained from that of Ox by 
an anticlockwise rotation through an angle a, Then the 
coordinates (x,y) of a point P with respect to the first 
system and its coordinates (x’,y') with respect to the 
second are related by 


z']_[ cosa sina] [x—a 
y'|~ |=sina cosa| | y—6]* 
Section 4: Square matrices and their inverses 


1. A square matrix is one with the same number of 
rows as columns. The main diagonal of a square matrix 
A consists of the elements a11,@22,... « 


2. A unit (or identity) matrix, written I, is a square 
matrix with elements 1 on the main diagonal and 0 
everywhere else. For any square matrix A and unit 
matrix I of the same size, we have 


AI=IA=A_ 


3. A singular square matrix is one whose rows are 
linearly dependent. A non-singular square matrix is one 
whose rows are linearly independent. Every non-singular 
square matrix A has a unique inverse A~' with the 
property 

AA SA wed 
Singular matrices (and non-square matrices) do not have 
inverses. 


4, To compute A~*, start with the matrix pair AjI and 
perform the same Gaussian elimination type row 
operations on both matrices, including multiplying or 
dividing a row by a non-zero constant. The linear row 
operations which convert the left-hand matrix A to a unit 
matrix I will convert the right-hand matrix I to A~'. 
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In particular, the 2 x 2 matrix 


as[r 


has inverse 
oe! d -b 
a -ae|4 ‘I. 


5. Matrix inverses have the property 
(Ar Aas: AntAay” = Ay ARS 


2 


6. A square matrix A is said to be symmetric if 
AT=A. 


Section 5: Introduction to determinants 


1, To every square matrix A there corresponds a number 
called its determinant, written det A or |Al. 


2. If Ais a 2 x 2 matrix, then 


ay a2 
a2) 022 


det A = 411422 — 012021. 


3. If A is a 3 x 3 matrix, then 
ay M12 a3 


detA =|a2 422 a25 
ay, 432 A393 
a22 «a23 @2,  aa3 
=a ~ a2 
432 ayy a3, ays 
42) a22 
tary : 
Qa) a32 


An n x n determinant can be defined in terms of 
(n —1) x (n — 1) determinants in a similar way. 


4, The rows (or columns) of a square matrix A are 
linearly dependent if and only if det A = 0. Therefore 
A~ exists if and only if det A # 0. 


5. Linear combinations of rows (or columns) of a 
determinant can be added to any other row (or column) 
without changing the value of the determinant. 


6. In an upper triangular matrix, all elements below 
the main diagonal are zero; in a lower triangular 
matrix, all elements above the main diagonal are zero; in 
a diagonal matrix al! elements not on the main 
diagonal are zero. The determinant of an upper 
triangular, lower triangular or diagonal matrix is the 
product of the elements on its main diagonal. 


7. To evaluate a determinant, use the operations of 
Gaussian elimination, namely adding or subtracting a 
multiple of one row from another and including row 
interchanges if necessary, to reduce the matrix to upper 
triangular form. Provided no row interchanges have been 
made, the upper triangular matrix has the same 
determinant as the original matrix; every interchange of a 
pair of rows reverses the sign of the determinant. 


Unit 21 Eigenvalues and eigenvectors 


Section 1: The theoretical eigenvalue problem 


1. Given a square matrix A and a number A, the 
equation 

Ax = Ax 
for the column vector x may have solutions other than 
the zero column vector. If it does, then we say that the 
number A is an eigenvalue of the matrix A, and that 
the non-zero column vector x is an eigenvector of A. 
Any non-zero scalar multiple of an eigenvector is also an 
eigenvector, so there are an infinite number of 
eigenvectors of A corresponding to each eigenvalue. 


2. One way of finding the eigenvalues A of a given square 
matrix A is to solve the characteristic equation 
det(A — AI) = 0, 
which is a polynomial equation in A. Eigenvectors can 
then be found by solving the equation 
(A-AI)x=0 
for x. 


3. If the characteristic equation has n different solutions 
Aayso+yAny where the size of the matrix A isn xn, then 
the eigenvalues of A are said to be distinct. If not, then 
the polynomial det(A — AI) has at least one repeated 
linear factor (A — \s)” with m > 2, and the 
corresponding eigenvalue 2, is said to be repeated, 


4, If the eigenvalues of an n x n matrix are distinct, then 
the matrix has n linearly independent eigenvectors. If the 
eigenvalues are not distinct, this may or may not be the 
case. 


Section 2: Iterative methods for finding 
selected eigenvalues 


1. If the square matrix A has eigenvalues 1, A2,---.An 

then for any real p,q: 

(i) A+qlI has eigenvalues A; +9,.-.,An +4: 

(ii) if A is non-singular, then AW" has eigenvalues 
Aitarestaet 

(ii) if A — pl is non-singular, then (A — pI)~* has 
eigenvalues (Ax ~ p)~*,.--,(An ~ p)~*5 

(iv) A, A+qlI, AW! and (A — pI)", if they exist, all 
have the same eigenvectors. 


2. Direct iteration. We use the recurrence relation 


Yrer = Ay,/Orsi, 
where a@,+1 is the element of largest modulus in Ay,. If 
Yo is not an eigenvector and not the zero vector, then as 
r becomes large, a, approximates to the eigenvalue of A 
with largest modulus (assumed real and distinct) and y, 
approximates to a corresponding eigenvector. 


3. Inverse iteration. We calculate y/.,, either from 
its definition y,,, = A~"y,, or by solving the 
simultaneous equations 


AY =¥rs 


and then obtain y,,, from 


Yost =Vr41/Ore1, 
where a,+1 is the element of largest modulus in y/,,,- 
For an arbitrary non-zero starting vector y, which is not 
an eigenvector, and with A non-singular, the number 
1/a, approximates to the eigenvalue with smallest 
modulus (assumed real and distinct) and y, approximates 
to a corresponding eigenvector, for sufficiently large r. 


4. Modified inverse iteration. To find the eigenvalue 
closest to p and a corresponding eigenvector of A, we 
‘start with an arbitrary non-zero vector yp which is not an 
eigenvector and calculate y’,,, either from its definition 
Yiu = (A —pl)~*y, or by solving 

(A= PDY 41 =¥or 
and then obtain y,,, from 


Yrs =Yrp1/Orsiy 
where a,+1 is the element of largest modulus in y/,,- 
If A — pl is non-singular, and r is large enough, the 
number p + (1/a,) approximates to the eigenvalue closest 
to p (assumed real and distinct), and y, approximates to 
a corresponding eigenvector. 


5. For a 3 x 3 matrix where we know the eigenvalues of 
largest and smallest modulus, denoted here by A; and Aa, 
a good choice of p in modified inverse iteration is either 
3(1A| + [Aal) oF —3(/Aa] + [Aal)- 


Section 3: Decomposition methods for finding 
all the eigenvalues 


1, To decompose a square matrix A is to express it as 
a product of two other square matrices, and the LU 
decomposition of A (when it exists) is 

A=LU, 
where L is lower triangular with all diagonal elements 
equal to 1, and U is upper triangular. 
In particular, a 2 x 2 matrix A has an LU decomposition 


ay aia] _ 1 0} }an a2 

aa, a22| [ar/ay 1) | 0 detA/an|* 
provided that a1, is not zero. 
2. The eigenvalues of a diagonal, upper triangular or 


lower triangular matrix are equal to the elements on the 
main diagonal. 


3. If A, B and P are n x n matrices, with P 

non-singular, then: 

(i) AB and BA have the same eigenvalues; 

(ii) A and P~'AP have the same eigenvalues. Further, 
if A has n distinct eigenvalues A1,.-.,An, and P is 
the matrix whose columns are the corresponding 
eigenvectors of A, then 


Ba 0 


PAP * 


Gi= == ss 
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where all off-diagonal elements of this last matrix 
are zero. Even if the eigenvalues are not all 
distinct, the result still holds provided a set of n 
linearly independent eigenvectors exists. 


4. The LR method for calculating eigenvalues. Given 
a square matrix A, form the sequence Ao, Ai,..., as 
follows: 

Ao=A, 

Ara = U-L,, 
where A, = L,U, is the LU decomposition of A,. For 
suitable A, the sequence Ao, Aj,... will converge to an 
upper triangular matrix whose diagonal elements are the 
eigenvalues of A arranged in order of decreasing modulus. 


5. Ifthe LR method applied to a matrix A fails to 
converge, try applying it to the matrix A + qI, where q is 
a real number. The eigenvalues of A + qI are 

Ar +q,A2 +4).++) where Ay, Ag,... are the eigenvalues 

of A. 


Unit 22 Simultaneous differential 
equations 


Section 1: First-order systems 


1, A linear first-order system of differential equations 
is one that can be written in the form 

Ai()x(t) + Ao(t)x(t) = g(t), 
where x(t) is a column vector whose elements are the 
unknown functions of t, x(t) is the column vector 
obtained by differentiating these functions, Ai(t) and 
‘A,(t) are given matrices and g(t) is a given column 
vector whose entries may depend on t. 
A linear second-order system is one that can be 
written 

A, (t)&(t) + Aa(t)x(t) + Aa(t)x(t) = g(t). 
If all the entries in the matrices A; and Ag, and As if 
present, are constants, we have a constant-coefficient 
system. If g(t) = 0 for all t, we have a homogeneous 
system. 


2. The normal form (if it exists) of a linear first-order 
constant-coefficient system is 

x(t) = Bx(t) + h(t), 
where B is a constant square matrix. 


3. If B is ann x n matrix with n linearly independent 
eigenvectors, say a1,-..,@n, corresponding to (not 
necessarily distinct) eigenvalues Ai,...,An, then the 
general solution of the homogeneous system 


x(t) = Bx(t) 


ME + Cnane’™’, 


where Ci,...,Cp are arbitrary constants. If B is real 
then the expressions a,e*"' in the general solution are 
either real or occur in complex conjugate pairs. If ae 
and @e™ is such a complex conjugate pair, then the real 
pair Re(ae*') and Im(ae**) may be used instead. 


x(t) = Ciae 
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Section 2: Further methods for linear 
first-order systems 


1. If B is an n x n matrix with n linearly independent 
eigenvectors, then the system 
x = Bx + h(t) 
can be solved by putting 
x= Py, 
where P is a square matrix whose columns are the 
eigenvectors of B. The resulting equation for y is 
y=P"'BPy +P™'h(t) 
and, since P~'BP is a diagonal matrix (its diagonal 
elements being the eigenvalues of B), we can solve the 
equations for the various elements of y individually. The 
solution x is then found from x = Py. 


2. Euler’s method, This numerical procedure can be 
used for first-order systems of the form 

x = m(x.t), 
with initial condition x(to) = xo, where m is a column 
vector whose elements may depend on the elements of x 
and on t. 


The approximation to x(t-), where t- = to +rh with h 
the step length, is denoted by 


ie aa ckaatee 

and we take Xo = xo. The recurrence relation is 
Xone = Xie bhm(Xryte) (F=1,....0)y 

or Xrp1 = Xp + hm(X,, tr). 


Section 3: Second-order homogeneous systems 


1. If B is an n x n real matrix with n linearly 
independent eigenvectors, the general solution of the 
system 


% = Bx 
is a sum of the following terms: 


(i) for each real positive eigenvalue A, with eigenvector 
ay, a term 


a,(CreV"" +. Dye V4); 
(ii) for a zero eigenvalue with eigenvector a,, a term 
ay(Cr + Det); 


(iii) for each real negative eigenvalue A, with 
eigenvector a,, a term 


a, (C;cos y/—Art + Dy sin y/—A,t); 
(iv) for each complex pair of eigenvalues A. and J, a 
term 


Cy Re(areY*"*) + Cl Im(are¥*"*) 
+D, Re(a,e~V*"") + Di Im(are~V*"), 
where a, is an eigenvector for Ar. 


Section 4: Forced oscillations 


1. Theorem 1: Suppose that x, is a particular solution 
of the system 


Ark + Ax + Asx = h(t). (1) 
Then the general solution of this system is 
X= Xp +X, 


where x. (the complementary function) is the general 
solution of the associated homogeneous system 


Aik + Aox + Asx = 0. 


2. Theorem 2: Suppose the matrices A1,A2,Ag are 
real, constant, symmetric and square, and that all their 
eigenvalues are positive. (A symmetric matrix is one that 
equals its own transpose.) Then every solution of 


Ai x(t) + Aox(t) + Asx(t) =0 
becomes small (i.e. every entry in the column vector x(t) 
becomes small) as t becomes large. 


3. To look for a particular solution of 
Aji % + A2x + Asx = bi coswt + basinwt, 
write the system in the phasor form 
Ark + Ark + Asx = Re(be™"), 
where b = b; — iba, and look for a solution of the form 
Re(ze"*"), 


where z is a column vector whose entries are complex 
constants. 


Unit 24 Normal modes 


Section 1: Modelling vibrating systems 


1. A lumped-parameter model is a model of a real 
mechanical system in which each component of the model 
has only one of the three properties mass, stiffness and 
damping. The only components used in this unit are 
particles (having mass only) and perfect springs (having 
stiffness only). 


2. The number of degrees of freedom in a 
lumped-parameter model is the number of coordinates 
necessary to specify the configuration (i.e. the positions 
of all the particles). 


3. The behaviour of a vibrating system is said to be 
forced if time-dependent external forces act on it; 
otherwise the behaviour is said to be free. 


Section 2: Free undamped vibrations with two 
degrees of freedom 


1. The equilibrium configuration (or static 
configuration) of a lumped-parameter system is the 
configuration for which the total force on each particle is 
zero; an equilibrium position of a particle is its 
position when the system is in an equilibrium 
configuration. 


2. A normal mode of a mechanical system is a motion 
of the system in which all the particles execute simple 
harmonic motion with the same angular frequency. The 
angular frequencies for which this is possible are called 
normal mode angular frequencies. 


3. In normal mode motion, the displacements of any pair 
of particles from their equilibrium positions are in a 
constant ratio, called the displacement ratio for that 
pair of particles. If the displacement ratio is positive, the 
two particles are said to be in phase; if it is negative, 
they are said to be phase opposed. 


4. The equations of motion for an undamped vibrating 
system can, by dividing each equation by the relevant 
mass, be put in the form 

% = Hx, 
where x is a column vector whose entries are the 
displacements of the various particles from their 
equilibrium positions, and H is a square matrix with 
constant entries. Then each eigenvalue of H is —w*, 
where w is a normal mode angular frequency, and the 
corresponding eigenvector, if the top entry in it is chosen 
to be 1, gives the displacement ratios for the 
corresponding normal mode. 


5. The most general motion of a vibrating system is an 
arbitrary linear combination of the normal mode motions. 


Section 3: Vibration absorption 


1. Forced vibrations. When a sinusoidally varying 
external force is applied to some part of a 
lumped-parameter system, a motion is in general possible 
for which all particles execute simple harmonic motion at 
the angular frequency of the external force. This 
steady-state motion can be used to model the 
steady-state response of a real system to such a force. 


2. If the frequency of the external force is close to a 
normal mode frequency, then the forced vibrations 
exhibit resonance—the amplitude of the steady-state 
response is very large. 


3. There may be a frequency for which one of the 
particles remains at rest in the steady state. This 
phenomenon can be used to isolate the object modelled 
by that particle from the effects of forced vibration at (or 
near) this frequency. The system is then said to act as a 
vibration absorber. 


Section 4: Extending the scope 


A degenerate system is one for which one of the normal 
mode angular frequencies is zero. In the corresponding 
motion each particle moves with the same constant 
velocity (rather than sinusoidally). 


Unit 25 Functions of more than one 
variable 


Section 1: First-order partial derivatives 


1. Any expression or formula which involves two 
variables, say x and y, and whose value is uniquely 
determined by the values of x and y, is called a function 
of the two variables x and y. If a variable = is a 
function of x and y (ie. : = f(x,y), where f is a function 
of the two variables x and y) then we call x and y the 
independent variables and = the dependent 
variable. The set of possible pairs of values (x,y) for the 
two independent variables is called the domain of the 
function; it can be thought of as a region in the 

(x, y)-plane. 

Functions of three (or more) variables are defined 
analogously. 


2. A function of two variables, say f(y), can be 
represented in three-dimensional space by the set of 
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points whose coordinates 2, y, = satisfy z= f(x,y). Such 
a set of points is an example of a surface. 


3. A section function is a function of one variable, 
obtained from a function of several variables by holding 
all but one of the independent variables fixed. 


4. If f(x,y) is a function of two variables, and a,b are 
two constants, then the (first-order) partial 
derivative of f(x,y) with respect to x at (a,b) is defined 
as the (ordinary) derivative at a of the section function 
f(x,b), and is denoted by Fa, b). Similarly, the partial 
derivative of f(x,y) with respect to y at (a,b) is the 
derivative at b (with respect to y) of the section function 


f(a.y), and is denoted by Hea, 6). 


5. We often use the partial derivatives of f(x,y) at (2,y) 
instead of at (a,b). The partial derivative Fria y) is 
calculated by differentiating f(x,y) with respect to x 
while treating y as a constant. This procedure is called 
differentiating f(x,y) partially with respect to x. 


Similarly, 5) is calculated by differentiating f(x,y) 


partially with respect to y at constant 2. 


First-order partial derivatives of functions of three or 
more variables are defined similarly. 


Section 2: Some applications of partial 
derivatives 


1, Let f(x,y) be a well-behaved function—that is, one 
such that the surface : = f(x,y) is smooth in the sense of 
having no folds or breaks. Then the Taylor polynomial 
of first order for f at (a,b) is 

of 


By (ar) 


miley) = Slab) +(e 0) 24 (0,8) + (yb) 


2. The first-order Taylor approximation formula at 
(a,b), f(z. y) ~ pi(x.y), holds for (x,y) close to (a,b). It 
can be written 


a 
Slew) ~ F(a,b) + (ea) 2L(a,6) + (y- eta, 6) 
oe ig it 
or fx oes Sty, 
where 
6f = f(x,y) = f(a,b), 
peas as 
eagle 


3. The triangle inequality is 
Ja +b] < [aj +b] 
and holds for any two numbers a and b. 


4. The tangent plane to the surface = = f(x,y) at the 
point (a,b, f(a,b)) has the equation = = pi(z,y)- 


5. The chain rule for the function of two variables 
f(z,u) is 

df _afdx , df dy 

dt Or dt” dy dt’ 
where x and y are functions of the single variable t. 


44 


6. The above formulas can be generalized to functions of 
three or more variables. 


Section 3: Higher-order partial derivatives 
and Taylor polynomials 


1. The second-order partial derivatives of a function 
f(z,y) are 


a7 ae (ae) 
dc? ~ Bz \Oz)" 
ia he aR (2) 
ayor ~ dy \ax)’ 
@f _ a (af 
Bray Ox (%) } 


ef a (#) 
Oy? Oy \ Ay ) © 
Third- and higher-order partial derivatives are defined 


similarly, as are the partial derivatives of functions of 
more than two variables. 


2. The Mixed Derivative Theorem. If f(x,y) isa 
well-behaved function of two variables, then 

Of  #F 

Oxdy — Aydx 


3. The Taylor polynomial of second order for 
Fa,y) at (a,b) is 


p2(z,y) = f(a,b) +(x — a) (a,6) +(y = 6) 54 (a,b) 


+3(2-a)5 FL (a, b) +(x —a)ly — )aray Zt (a,b) 


+ sy Lab). 


4. For (x,y) near to (a,b), the second-order Taylor 
approximation at (a,b) is f(«,y) ~ pa(x,y) and is 
normally more accurate than the first-order Taylor 
approximation. 


Section 4: Maxima and minima 


1. If, for all x and y sufficiently close to a and b, 


S(a,y) < f(a,), 
then f has a (local) maximum at the point (a,b); but if 
the inequality is 

f(x,y) = F(a,b), : 
then f has a (local) minimum at (a,}). 


2, An extremum is a point which is either a local 
maximum or a local minimum. 


3. A stationary point of a function f(x,y) is a point 
(a,b) where 


Fate Lipase 
Be) =0 and 5 (a,b) =0. 


‘The stationary point criterion. If (a,b) is an 
extremum of a well-behaved function f(x,y), and is not 
on the boundary of the domain of f(z, y), then (a,b) is a 
stationary point of f(z,y). (The converse statement is 
not necessarily true.) 


4, A saddle point is a stationary point which is not an 
extremum. 


5. The AC — B? criterion. Suppose (a,b) is a 
stationary point of f(x,y). oH let 


a= Loan, a= gf (a,b), c= Shea »). 


Then 
maximum 
minimum 


if A <0, 


(i) it AC ~ BY > 0, (0,8) ina { if A > 0; 


f AC — B® <0, (a,b) is a saddle point; 
if AC — B? 


(ii) 


(iii) = 0, the criterion gives no information. 


Unit 26 Vector calculus 


Section 1: Scalar and vector fields 


1. A scalar field function, or scalar field, is a 
function which to each point in a given region of space 
associates a unique scalar. The region is the domain of 
the scalar field. 


2. The contour curves (or contours) of a 
two-dimensional scalar field @(z, y) are the curves 
(x,y) = constant. The contour surfaces of a 
three-dimensional scalar field $(z, y, 2) are the surfaces 
(x, y,2) = constant, 


3. A vector field function, or vector field, is a 
function which to every point in a given region of space 
(or of the plane) associates a unique vector. 


4. The field lines of a vector field are the family of 
continuous curves such that at every point the tangent to 
the curve is along the direction of the vector field at that 
point. 


Section 2: Differentiating scalar fields 


1. The symbol V, jhe by 
a 
v= ig; +35 + KZ, 
ix 


isa differential vector operator. 


2. For a three-dimensional scalar field ¢(x,y, =), the 


is called the gradient of ¢ and is denoted by grad @ or 
Vo. 


3. The direction of grad é at a point P in the domain of 
the scalar field @ is perpendicular to the contour curve or 
surface passing through P. This direction is the normal 
to this contour at P, 


4. The vector field grad @ gives the magnitude and 
direction of the maximum (spatial) rate of change of the 
scalar field @ at any point P in the domain of the scalar 
field. 


5. The (spatial) rate of change of a scalar field ¢ at a 
point P in the direction of a unit vector e is given hy 
e- (grad @). 


Section 3: The scalar line integral 


1. The scalar line integral of a vector field F along a 
curve C, called the path of integration, with 
end-points A and B, is defined as 


F dr = lim) F, -6r,, 
eo oa ae 


= 
where the curve is divided into n segments with 
end-points P,(= A), P2, Ps,..-.Pn+i(= B); F; means 


F(P,); 6r, means P;P;+1; and the limit is taken in such a 
way that n becomes very large and the segments all 
become very short. The end-points A and B may 
coincide, in which case the path of integration is closed; 
if they do not coincide, it is open. 


2. A three-dimensional curve can be expressed in terms 
of a parameter ¢ by the parametric equations 
w=a(t), y=y(t), z=2(t) (@st<b). 


3. To evaluate a scalar line integral [ F-dr: 
ic 


(i) describe the path C by means of a set of parametric 
equations; 
(ii) use these parametric equations and the formula 
dr _ de, dy, de 
a at ata 
to express the vector field F and dr/dt in terms of 
the parameter t, and hence rewrite the line integral 
as an integral with respect to t: 


F-de= [ F(t). “at; 
Pessoal fe ae 


(iii) evaluate this integral with respect to t. 


4. The line integral of grad @ is independent of the path 
apart from its end-points: 


(grad @) - dr = 6(B) — (A). 
AB 


5. A vector field F is conservative if for every closed 
path C in the domain of F we have 


fF-ae=0, 
ic 


where the symbol f 


Cc 
closed path C. Otherwise the vector field is called 
non-conservative. 


means a line integral round the 


Section 4: The curl of a vector field 


1. For a vector field F we define the component of 
curl F at a point P in the direction of a unit vector n to 
be 

1 


a-(curlF) = fim 3 fdr, 
aoe Cc 


where C js a small closed curve in the plane through P 
perpendicular to n, which contains P and encloses an 
area A. (The direction of integration round C is chosen so 
that a screw turned in this direction advances along n.) 


2, The curl of a vector field F(x, y,2) = Fii+ Pej + Rik 
can be evaluated by using the derivative formula 


dy Oz 


OF: _ oF, 
+(2 as a) 


culF=VxF= (B-F) 


3. The curl of a vector field F is a detector for rotational 
motion; for instance, the curl of a fluid velocity field 
describes the spin of a small float placed in the fluid. 


4. If F is a conservative vector field: 

(i) f F - dr = 0 for every closed path C in the domain 
of F; 

(ii) the line integral F - dr is independent of the 


AB 
path taken between any two points A and B; 
curl F = 0; 

(iv) there exists a scalar field 6 such that F = grad ¢; 


(v) / F - dr = (B) — ¢(A), where F = grad. 
AB 


5. Provided that the domain of F is sufficiently simple in 
shape, the converse of property 4(iii) above is also true, 
ie. the equation curl F = 0 implies that F is conservative. 


Unit 27 Multiple integrals 
Section 1: The surface integral 


1. Suppose that f is a two-dimensional scalar field whose 
domain includes a region S of the (2, y)-plane. We divide 
the region S into a number N of area elements so that 
the ith element has area 5A;, and denote some point in 
the ith element by P;. Then the surface integral of f 
over the region S is defined to be 


* 
[444 = tim 107 )64. 
s 1 


where the limit is taken in such a way that N becomes 
very large and the sizes of all the elements become very 
small, We call S the region of integration. 


2. The following procedure can be used to evaluate 

surface integrals. 

(i) Draw a diagram showing the region of integration, 
including the equations of its boundary lines. 


(ii) Mark on the diagram the minimum and maximum 
values of x, which we call here a and 6 respectively. 


(iii) Draw a typical strip of the region parallel to the 
y-axis and denote the formulas for the minimum 
and maximum values of y (which depend on the 
x-coordinate of the strip) by a(x) and (x) 
respectively. 

(iv) The surface integral can now be written in terms of 
two single integrals: 


frees fo [tone 


(v) Evaluate the ‘inner’ integral, holding x constant: 


(vi) Evaluate the ‘outer’ integral: 


=b 
[sa-f g(x) dx. 
Ss p= 


An analogous procedure, using strips parallel to the 
z-axis, can also be used in which we first integrate 
over x and then integrate over y. 


Section 2: Multiple integrals useful in 
mechanics 


1. Suppose that f is a three-dimensional scalar field 
whose domain includes a region B of space. The volume 
integral of f over the region B is defined as 


N 
i fadV =lim > f(P)Vi, 
B 1 


where B is divided into a number N of volume elements, 
the ith of which has volume 6V, and contains the point 
P;, and the limit is taken in such a way that all the 
volume elements become very small. We call B the 
region of integration. 


2. If pis the scalar field giving the local density of a 
body occupying a region B, then the mass of the body is 


M [oa 
B 


The moment of inertia of this same body about some 
given axis is 


r= [ (otyav, 


where d denotes distance from this axis. 


3. The mass of a hollow sphere of uniform density p with 
external radius a and internal radius b is 4xp(a° — b°), 
and its moment of inertia about an axis through its 
centre is Amp(a° — 6°). The corresponding formulas for a 
solid sphere are obtained by setting b = 0. 


4. To evaluate a surface integral if dA using polar 
s 
coordinates (r,8) the following procedure can be used. 


(i) Draw a diagram showing the region of integration, 
S, and its bounding lines. 


(ii) Write the surface integral in terms of r and @ using 


the formulas 
=rcosé, 
y=rsin@, 
6A =r br 60. 


In the limit of very small area elements, this last 
formula becomes dA = r dr dé. 


(iii) Find the minimum and maximum values of @ in S, 
say 8; and 62. 

(iv) For a strip of fixed 8, indicate the minimum and 
maximum values of r, say a and b (here assumed 
independent of , for simplicity), 


(v) Write the surface integral in terms of two single 
integrals: 


8=82 [ pr=b 
dL [ F(r.8) rar| 8. 
o=0, rae 


(vi) Evaluate the inner integral, then the outer integral. 


Section 3: The volume integral 
The following procedure can be used to evaluate a 


volume integral | f dV using Cartesian coordinates, 
B 

(i) Draw two diagrams, showing (a) the region B of 
integration and (b) the projection S$ of this region on 
the (x, y)-plane. (The projection of a point (x,y, =) 
on the (x,y)-plane is the point (x,y), and the 
projection of B comprises the projections of all the 
points in B.) 

(ii) Within the region B, draw a column parallel to the 
z-axis and mark its end-points, which are the 
equations of the top and bottom surfaces of B 
written in the form z= (x,y) and z = o(z,y) 
respectively. 


(iii) Evaluate the single integral 


e=w(Ey) 
g(x,y) -[ S(x,y,2) dz. 


=O ey) 
(iv) Evaluate the surface integral of g(z, y) over the plane 


region S, The result is nf fav. 
a 


Similar procedures, with different orders of integration, 
can also be used. 


Section 4: Changing variables in three 
dimensions 


1. The cylindrical polar coordinates (r,4, =) of a 
point P in space are related to its Cartesian coordinates 
(z,y, =) by 


z=rcosé, 
v=rsin@, 
2=2, 


where r > 0,0 <0 < 2m, —c0 < = < 00. That is, (r,0) are 
the plane polar coordinates of the projection of P onto 
the (x,y)-plane. 


The volume element in cylindrical polar coordinates is 
8V =r br 6062. 
This becomes dV = rdr d@ dz in the limit. 


2. To evaluate a volume integral using cylindrical polar 
coordinates, a procedure similar to the one given above 
for Cartesian coordinates can be used; the only difference 


is that we use plane polar coordinates r,@ in place of 
plane Cartesian coordinates x,y throughout. 


3. The spherical polar coordinates (r,6,¢) of a point 
P are related to its Cartesian coordinates (x,y,z) by 

x =rsinOcos@, 

y=rsin@sind, 

z=rcos6, 
where r > 0,0 < @ < 24,0 <0 <7. In this case r is the 
distance from P to the origin O, @ is the angle between 
the line OP and the positive z-axis, and ¢ is the angle 
between OQ and the positive x-axis, where Q is the 
projection of P onto the (z, y)-plane. 


cj Q 
The volume element in spherical polar coordinates is 
6V =r? sin O br 6060. 


This becomes dV = r? sin @ dr d0 d@ in the limit. 


4. To evaluate a volume integral using spherical polar 
coordinates, use a procedure similar to the one given for 
Cartesian coordinates by expressing the volume integral 
as three repeated single integrals. 


Unit 28 Moments and circular motion 


Section 1: The components of a force 


1. The component of a force F in a direction making an 
angle a with that of the force is equal to F cosa, where 
F=|F|. 


2. Ifa particle is static (in equilibrium), then the sum of 
the components of the forces acting on the particle in any 
direction is equal to zero. 


Section 2: The moment of a force 


1. The magnitude P of the moment of a force about a 
fixed point O is the product of the magnitude F of the 
force and the perpendicular distance d from the point O 
to the line of action of the force, that is, 


lr=Fd. 


2. The total moment of the forces acting in two 
dimensions on a body about any point is the sum of the 
anticlockwise moments minus the sum of the clockwise 
moments, 


3. Ifa rigid body is static (in equilibrium), then both the 
total force on the body and the total moment of forces 
acting on the body about any point are zero. 


aT 


Section 3: Motion in a circle 


1. In terms of polar coordinates, the radial and 
transverse components of the position, velocity and 
acceleration of a particle moving in a circle, r = constant, 
are as follows. 


Radial | Transverse 

component | component 
Position r 0 
Velocity 0 a) 
Acceleration | —ré” ro 


2. For circular motion with radius r and angular 
velocity w = @, the velocity of the particle is v = rw, 
tangential to the circle. Both v and w are positive for 
anticlockwise motion. The component of acceleration 
towards the centre of the circle is rw? = v?/r. For 
non-uniform motion, the acceleration also has a 
transverse component rw. 


3. For uniform circular motion, with radius r and 
constant angular velocity w, the acceleration has 
magnitude ru* = v?/r and is directed towards the centre 
of the circle. The time taken for one complete revolution 
of the circle is tr = 2n/|w]. 


Section 4: Using vector notation 


1. The torque I, about a fixed point O, of a force F 
acting at a point which has position vector r relative to 
the point O, is defined to be 


T=rxF. 


2. The torque of all the gravitational forces acting on a 
system of particles or extended body can be evaluated by 
considering the total weight of the system to be 
concentrated at its centre of mass. 


3. If a system of particles or extended body is in 
equilibrium, then the total external force is equal to zero, 
as is the total torque of the external forces about any 
fixed point. 


4. With respect to plane polar coordinates [r, 6], the unit 
vectors e, and eg (in the radial and transverse directions, 
respectively) are given by 

e, =cos@i+sin0j, e9 =—sin#i+cos0j. 
The vectors e,,es,k form a right-handed triad of unit 
vectors. The directions of e, and eg vary with time, so 
that 


@,=6es and 


5. The position, velocity and acceleration of a particle 
moving in a circle of constant radius r in the (x,y)-plane, 
with centre at the origin, are 


r= 


er, 


v 


and a=i=—ré’e, + rées. 
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6. The angular velocity w for a particle in circular 
motion is a vector whose magnitude is equal to the 
angular speed, and whose direction is along the axis of 
rotation, in the sense in which the rotation would drive a 
right-handed screw. If the unit vector k is perpendicular 
to the plane of motion and the centre of the circle lies on 
the z-axis, then 


w= dk. 


7. For circular motion with angular velocity w, the 
velocity of a particle is given by 

r=wxr, 
where r is the position vector of the particle relative to 
an origin on the axis of rotation. 


Unit 29 Angular momentum and rigid 
bodies 


Section 1: Angular momentum 


1. The strong form of Newton’s third law states 
that for a pair of interacting particles, the force Fy2 on 
Particle 1 due to Particle 2 is equal in magnitude but 
opposite in direction to the force Fz; on Particle 2 due to 
Particle 1, that is, 


Fi2 = -Fa, 


and each of these forces acts along the line which joins 
the two particles. 


2. For a particle which has position vector r relative to 
an origin O and linear momentum p = mf, the angular 
momentum [ relative to O is 


l=rxp=rxmi. 


3. The torque law for a single particle states that the 
rate of change of the particle's angular momentum 
relative to a fixed point is equal to the applied torque 
relative to that point, that is, 

i=4. 
In the special case where the total torque acting on the 
particle is zero, the angular momentum of the particle 
about the fixed point is constant, In other words, if y= 0 
then lis constant. This is the law of conservation of 
angular momentum for a single particle. 


4. The total angular momentum L of a system of n 
particles, relative to a given origin, is the vector sum of 
the angular momenta of the individual particles, that is, 


L= Sau = ae x MAFi), 
=I =I 


where r; is the position vector of the ith particle and m, 
is its mass. 


5. The total external torque P'* acting on a system 
of n particles, relative to a given origin, is defined to be 


the vector sum of the external torques which act upon 
the individual particles, that is, 


rn Saft 


6. The torque law for a system of particles states that 
the rate of change of the total angular momentum about 
the origin is equal to the total external torque about the 
origin, that is, 

Lao 
This law applies to any system for which the inter-particle 
forces satisfy the strong form of Newton's third law. 
In particular, if the total external torque is zero then the 
total angular momentum L is constant. This is the law 
of conservation of angular momentum for a system 
of particles. 


Section 2: Rigid-body rotation 


1. A rigid body is a many-particle system with the 
property that all the inter-particle distances remain 
constant in time. The inter-particle forces in a rigid body 
are assumed to satisfy the strong form of Newton's third 
law, so that the torque law L = P*** applies. 


2. Rigid-body rotation about a fixed axis is a 
rigid-body motion for which particles on the axis remain 
fixed, while other particles move in circles centred on and 
perpendicular to the axis. If the axis of rotation is chosen 
to be the z-axis, then all particles have the same angular 
velocity w = 6k. If r; is the position of Particle i, then its 
velocity is given by 


Ok & ry. 


R=wxrn= 


3. The 2-component of the total angular momentum L of 
the body described above is 


L, = Ib, 


where 


t= Siaut 
is 


is the moment of inertia of the body about the axis of 
rotation (the sum being taken over all particles of the 
body). It follows from the torque law that 


lo=L, =r. 
‘The kinetic energy of the body is given by 
T= 16". 


4. The moment of inertia of a continuous rigid body 
about a given axis of rotation is 


= [ p(r)d?(r) dV, 
B 


where p(r) is the density of mass at point r, d(r) is the 
distance of the point r from the axis of rotation, and the 
integration is over the volume of the body B. 


For a homogeneous body, the density p is constant. The 
moments of inertia of certain homogeneous rigid bodies 
about axes through their centres of mass are given in 
Table 1. 


Table 1: Moments of inertia of homogeneous rigid bodies 


Homogeneous rigid Diagram Dimensions | Moment of inertia 
body of mass M; about given axis 
axis 

solid sphere; radius R 2MR? 


axis through centre 
of mass 


hollow sphere; 
axis through centre 
of mass 


inner radius a, 
outer radius R 


iu (E=S) 


~2MR if R~a 


solid cylinder; 

axis through centre 
of mass and along 
axis of cylinder 


radius 


i MR? 


hollow cylinder; 
axis through centre 
of mass and along 
axis of cylinder 


D|( | |e) 


: 


inner radius a, 
outer radius R 


4M(R? +a?) 
= MR ifR~a 


solid cylinder; 
axis through centre 
of mass and 
perpendicular to 
axis of cylinder 


radius R, 
length h 


IMR? + Mh? 
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Table 1 (continued) 


perpendicular to 
one pair of faces 


ae 
Loa 


— 


Homogeneous rigid Diagram Dimensions | | Morent ofimertn 
body of mass M; about given axis 
axis 
hollow cylinder; inner radius a, | +M(R? + a?) + 4Mh? 
axis through centre outer radius R, 
of mass and length h ~ $MR? + 4Mh? 
perpendicular to R ifR~a 
axis of cylinder 

a 
solid rectangular | faces EM(a? +84) 


brick; perpendicular 
axis through centre to axis have 
of mass and sides of 


lengths a and b 


thin rod 
(cross-section of 
arbitrary shape); 
axis through centre 
of mass and 
perpendicular to 
rod 


h 
Rt 


a = 


length h 


ayo 


2 


5. The parallel axes theorem states that the moment 
of inertia of a rigid body about an axis A can be written 
as 

1=Ia+Md’, 
where M is the mass of the body, d is the distance of the 
centre of mass G from the axis A, and J¢ is the moment 
of inertia of the body about an axis B which passes 
through G and is parallel to A. 
Section 3: Rotation about a moving axis with 
fixed orientation 


1. Consider a system of n particles, with centre of mass 
R, total mass M, total kinetic energy T and total angular 
momentum L about the origin, which is acted upon by a 
total external force F* corresponding to a total external 
torque P** about the origin. Suppose that the ith 
particle has mass mj and position r;, and experiences an 
external force F¢**. Then the position of the particle 
relative to the centre of mass, r**!, is defined by 

lon -R. 
The centre of mass decomposition theorems are as 
follows. 


(i) The total kinetic energy of the system is given by 
T =To+}M\|R/, 


where = 
arely2 
To = £ > miles) 
‘=n 
is the kinetic energy relative to the centre of mass. 


(ii) The total angular momentum of the system about 
the origin is given by 
L=Lo+Rx MR, 
——where 


Le = Sor" x mia 
rast 


is the angular momentum about the centre of mass. 
(iii) The total external torque of the system about the 
origin is given by 
rt =fe+RxF™, 
where 


Po= or x Fe 
‘=I 


is the total external torque about the centre of 
mass. 


2. For a rigid body rotating with angular velocity w = 0k 
about an axis which passes through the centre of mass 
and is always aligned with the z-axis, the kinetic energy 


and the z-component of the total angular momentum 
relative to the centre of mass are given respectively by 


Te = 310° and Le-k= 10, 


where J is the moment of inertia of the body about the 
axis of rotation. 


3. For any system of particles, the rate of change of the 
total angular momentum relative to the centre of mass is 
equal to the total external torque relative to the centre of 
mass, that is, 

Lc =a. 
This is the torque law relative to the centre of 
mass, 


4. Ifa cylinder or sphere of radius R rolls (without 
slipping) along a rough plane, then the speed w of the 
centre is related to the angular speed |6| of the body 
about its centre by the rolling condition 


v = Rij. 


Section 4: Rotation about an axis whose 
orientation changes 


1. Ifa body spins around an axis of symmetry (the spin 
aris), and this axis is fixed, then the total angular 
momentum of the body is given by 

L=J6n, 
where / is the moment of inertia with respect to the spin 
axis and @n is the angular velocity of the body about this 
axis. 


2. Ifa body spins around an axis of symmetry which is 
not fixed, and J, én are defined with respect to the spin 
axis as above, then the approximation 

L~Jén 
is valid provided that |1| is small compared with 0. 


3. The angular momentum of a spinning top has a 
constant magnitude and a constant vertical component. 
It satisfies the equation 
L=akxL, 
where & is the rate of precession of the spin axis about 
the vertical axis, and k is a unit vector in the vertically 
upward direction. If the top has mass M, moment of 
inertia J and angular velocity én about the spin axis, 
then the rate of precession is given by 
Mage 
a>, 
10 
where c is the distance of the centre of mass from the 
point of contact with the ground. 


Unit 30 Planetary orbits 


Section 1: Review of energy and planar 
motion 


1. For particle motion in three or fewer dimensions, a 
force F acting on the particle is conservative if there 
exists a scalar field function U = U(r) = U(z,y,z) such 
that 


—gradU. 
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The function U is the potential energy of the particle. 
Tn such a case, the total mechanical energy of the particle, 


E= 
(where #? = #-# = ||*), is constant throughout the 
motion. 


2. In one dimension, the condition for a conservative 
force and the expression for the total mechanical energy 
become respectively 
dl 3 
=- and B= mz? +U, 


where F is the z-component of the force and U = U(z). 


3. The unit vectors e,, eg corresponding to the plane 
polar coordinates |r, @] (where x = rcos@, y = rsin®) are 
related to the Cartesian unit vectors i,j by the equations 


e, =cosBi+sin8j, e9 = —sin0i+cos8j 
and i=cos#e,—sin@es, j=sinde, +cosbey. 


4. The derivatives é-,@9 are given in terms of e,,e9 by 
@,=Bey and & =—ée,. 
5. The position r, velocity # and acceleration # of a 


particle moving in the (,y)-plane are given in terms of 
plane polar coordinates by the equations 


r=rer, 
F=fe, +réeg, 


#=(7-ré")e, + 14 (ré)eu. 
Section 2: Kepler’s laws and Newton’s law of 
gravitation 


1, Kepler’s laws of planetary motion are as follows. 


(i) Each planet moves in an ellipse, with the Sun at 
one focus, 

(ii) The line joining a planet to the Sun sweeps out 
equal areas in equal times. 

(iii) The square of the orbital period of a planet is 
proportional to the cube of the semi-major axis of 
its orbit. 


2. (i) The ellipse shown in the figure has semi-major 
aris a and semi-minor aris b. With respect to the 
axes O'XY, its equation is 

Mee Ae 
ate 
where a > 6b > 0. 


=1, 
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(ii) ‘The area of the ellipse is rab. 


(iii) The foci of the ellipse are located a distance q to 
either side of the centre O’, where 


q= Vat 
(iv) With respect to axes Oxy with origin at the 


right-hand focus, and in polar coordinates, the 
equation of the ellipse is 


1 1 4 ec0s8, 
: 
where 
2 
tals 
& 


is the semi-latus rectum, and 


q [_e 
=> aa 


is the eccentricity, where 0 < e < 1. 


3. Newton’s universal law of gravitation states that 
the gravitational force on a particle with mass m, and 
position vector r; due to another particle with mass m, 
and position vector ry is 


where G is the gravitational constant. 


4. A body has a spherically symmetric mass 
distribution if (with the origin of coordinates at the 
centre of the body) its mass density function p(r) 
depends only on the radial distance 

|r| =r = \/x? + y? + 27. The gravitational force exerted 
by such a body B on an exterior particle with mass m 
and position vector r is given by 


Win} Seat 


7 
where e, = r/r is the unit vector in the outward radial 
direction, 


is the gravitational potential energy of the particle, and 
M= [i pir’) dV 


B 
is the total mass of the body. 


Section 3: Orbits 


1. A central force F, acting from the origin on a 
particle with position vector r, has the form 
F(r) = g(r)r, 
where g(r) = g(x,y, 2) is a scalar function of position. If 
the particle has mass m and is not acted upon by any 
other force, then it has constant angular momentum 
L=rx mi. If the unit vector k is chosen in the 
direction of L, then 
(i) the motion of the particle is confined to the 
(x,y)-plane; 
(ii) L= Lek, where L. = mr?@ is constant. The 
gravitational force due to a spherically symmetric 
body whose centre is fixed at the origin is a central 


force. In this case, the constancy of L. is equivalent 
to Kepler’s second law. 


2, An isotropic central force has the form 
F(r) = f(r)er, 


where r = \/x? + y? + 2? is the distance from the origin, 
and e, =r/r. For a particle of mass m acted upon by 
such a force, the radial equation of motion is 
- Le 

mr = f(r) + ane 
which may be integrated (analytically or numerically) to 
give r(t). The gravitational force is isotropic, with 
f(r) =—GmM/r?, where M is the central mass and G is 
the gravitational constant. (It is assumed that M is much 
larger than m, in order that the central mass should 
remain fixed at the origin.) 


3. For the particle moving under gravity as described 
above, the total mechanical energy E is constant. It may 
be expressed as 
(i) E = tm(i? +H) + U(r), 
where U(r) = —GmM/r is the (actual) potential 
energy of the particle; 


(ii) E= im +u%(r), 
where 
pie eechEs Gm 
WO = seed 


is the effective potential energy. Orbits are 
unbound if £ > 0, and bound if Ug, < E <0, 
where USM, is the minimum value of U(r). 


4. The equation of gravitational orbit for the particle 
described above is 


Fai +ecose, 
where the semi-latus rectum / is given by 
B 
'= MG" 


If 0 < ¢ < 1 then the orbit is an ellipse with eccentricity 
¢, showing that Kepler's first law is satisfied. If ¢ > 1 
then the motion is unbound. 
In terms of the eccentricity and the semi-latus rectum, 
the total mechanical energy is 

GmM 

E= 
2 


(e? —1). 


5. Kepler's third law is also a consequence of this 
Newtonian model, since the period T of an elliptical orbit 
and its semi-major axis a are related by the equation 
4n* 
T? =ka*, where k= —_. 
ka*, where k G 

Section 4: Central forces and conservation 
laws 


1, Ifa scalar function ¢(r) depends upon position only 
through the radial distance |r| =r = (x? + y? + 2*)'/?, so 
that (r) = @(r), then the gradient of ¢ can be written as 


where e, = r/r is the unit vector in the outward radial 
direction. 


2. Hence if F is an isotropic central force, that is, 
F(r) = f(r)er 


for some function f, then F is conservative, with a 
potential energy function U(r) such that 
a 
f(r) =- ar 
The motion of a particle in this force field is planar, and 
has constant total mechanical energy 
B= 3m(# +76") + U(r). 


Unit 31 Fourier analysis 


Section 1: What Fourier analysis is used for 


1. A function f(t) is said to be periodic, with period 
T, if 

F(t)=f(t+nT) for allt (n= +1,+2,+3,.. 
that is, the graph of f(t) consists of repeating blocks of 
length T, We usually take T as the size of the smallest 
such block, but any multiple of this also satisfies the 
definition of a period. 


2. If the response of a mechanical system to an input 
Fosinwt is Asin(wt +), we call A/Fo the relative 
response and ¢ the phase shift. The frequency 
response of the system is described by giving A and ¢ as 
functions of w. For a linear system (i.e. one that can be 
modelled by a linear differential equation whose input is 
the right-hand side of the equation and whose output is 
the unknown function), the principle of superposition (see 
Unit 6) enables us to find the output produced by any 
input that is a linear combination of sinusoidal functions. 
That is, if C,,C2,... are constants and F,(t), F2(t),... 
are sinusoidal functions, and if the input F,(t) acting 


alone produces output z-(t), then the input }~> C-F,(t) 
Fi 
ral 


produces the output } > Cyar-(t). 


3. The sum of a number of sinusoidal terms of different 
frequencies is a non-sinusoidal function and in general is 
also non-periodic; but if the frequencies of the sinusoids 
are all multiples of w, then their sum will be periodic and 
will have period 2x/w. In this case we call w the 
fundamental angular frequency. 


Section 2: Fourier series for periodic functions 


1. If f(t) is a periodic function with period T, then the 
Fourier series for f(t) has the form 


2 
F(t) = M+ )\(Ansinnwt + By cos nut), 
oat 


where w = 2n/T is the fundamental angular frequency. 
The constants in this series are given by 


1 st? 
M=—> t) dt, 
z/ we 


1, 


2 7? 
an= 3 f f(t)sinnwtdt (n =1,2,3,...), 
“7/2 
2 pt 
a= 2 f F(t)cosnutdt (n= 1,2,3 
-r/2 


2. A function is said to be continuous if its graph is an 
unbroken curve. A function is piecewise continuous if 
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its graph consists of continuous pieces, with only a finite 
number of such pieces on any finite part of the domain 
and each piece having a well-defined end-point. 


3. Theorem 2.3: If a periodic function and its 
derivative are piecewise continuous, then the function is 
equal to its Fourier series for all values of t at which there 
is no gap in the graph. 


Section 3: Some outstanding points 


1. A function f(t) is said to be even if 
F(t) = f(—t) (for all t). 
The graph of an even function is symmetric under 


reflection in the vertical axis. The Fourier series for an 
even function contains no sine terms. 


An odd function f(t) is one satisfying 
f(t) =—f(—-t) (for all t). 
The graph of an odd function is symmetric under a 


rotation about the origin through 180°, The Fourier series 
for an odd function consists of sine terms only. 


2. We can often obtain useful approximations to a given 
periodic function f(t) by truncating its Fourier series, i.e. 
using a limited number of terms. 


Section 4: Extending the scope 


1. The theory described in this unit is not restricted to 
periodic functions of the variable t; any variable, for 
example x, may be used. 


2. A non-periodic function whose domain is a finite 
interval of the t-axis can be represented by a Fourier 
series; we do this by extending the definition of the 
function to the entire t-axis in such a way that the 
extended function is periodic, and then computing the 
Fourier series for the extended function. This series then 
represents the original function on its own restricted 
domain. 

If the Fourier series is required to have sine (or cosine) 
terms only, it is necessary to use a periodic extension 
which is odd (or even). 


Unit 32 Partial differential equations 
Section 1: Introduction and definitions 


1. A partial differential equation is an equation 
involving one or more of the partial derivatives of some 
unknown function of two or more variables, and possibly 
also the value of the function itself and the values of the 
independent variables. (In contrast, an ordinary 
differential equation is a differential equation where 
the unknown is a function of just one variable, and which 
therefore does not involve partial derivatives.) Solutions 
of partial differential equations may be required to satisfy 
boundary conditions on the boundary of the domain 
of the unknown function. If one of the independent 
variables represents time, a condition imposed at the 
initial time is usually called an initial condition rather 
than a boundary condition. 
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2. An anes of the form 


aes oa + it pees jh 


+d(zx, Be e(x, 2 + f(z,t)U = 9(x,t), 


Se Oz 
where U(z, t) is the erly function and at least one of 
the given functions a,b,c is not the zero function, is 
called a linear second-order partial differential 
equation. If in addition g(zx,t) is the zero function, we 
say the equation is homogeneous; if the functions 
a,b,...,f are all constants, we say the equation has 
constant coefficients. 


3. Theorem 1: If u and v are two solutions of any 
homogeneous linear partial differential equation, and A 
and B are arbitrary constants, then Au + Bv is also a 
solution. 


Section 2: Properties of coiled springs and 
elastic strings 


1. A spring is in a state of tension if it is extended, 
and in a state of compression if it is compressed. The 
modulus of elasticity k of a spring is defined by 

K 

hy’ 

where K is its stiffness (denoted by k in Unit 7) and lo is 
its natural length. If a spring, or string, is in a state of 
tension, then for every point P along its length we define 
the tension T at P to be the magnitude of the force 
exerted by the part of the spring to one side of P on the 
part to the other side of P. The value of T is always 
positive for a spring or string in a state of tension. 


ea 


2. For a spring stretched between two fixed points A, B, 
we define the reference state as the configuration which 
makes the tension the same all along the spring. For a 
general state of the spring, still confined to the line AB, 
the tension T at a pia P on the spring is given by 


= (To+ nat —h, 
where 
(i) Tp is the tension in the reference state, 
(ii) kis the modulus of elasticity, 
(iii) a is the position of P in the reference state, 
measured with respect to an x-axis lying along AB, 


(iv) s(x) is the position of P in the general state. 


If the spring is moving, so that s depends on the time ¢ 


as well as on 2, then e in the above formula is replaced 


a 
by Feltt)- 


3. A spring whose mass is not negligible is called a 
heavy spring. Vibrations of a stretched heavy spring in 
which all motion is along the direction of the spring itself 
are called longitudinal vibrations; vibrations in which 
the motion is at right angles to this direction are called 
transverse vibrations. Longitudinal vibrations can be 
modelled by the equation 

eu 1 eu 

02 2 OP 
where 


(1) 


@) U=s-—z, 
(ii) = (To+k)/m, 


(iii) _m is the mass of the spring per unit length in the 
reference state. 


If the ends of the spring are fixed at the points 
corresponding to z = 0 and x = L, where L is given, the 
boundary conditions are 
U(0,t) =0 
UL) =0 f Orallt. (2) 


Equation (1) is called the wave equation and the 
constant ¢ in it is called the wave speed. 


Section 3: The method of separation of 
variables 


1. The method of separation of variables, applied to 
the wave equation (1) with boundary conditions (2), is to 
look for solutions of the form 

U(x,t) = X(x)T(t), (3) 
where X,T are functions (of one variable) to be found. 
Substituting into the wave equation (1) and dividing by 
X(zx)T(t) gives 

Xo) _ Veet 

X@) FT) — 
where 1 is a constant. So the function X satisfies the 
ordinary differential equation 


(4) 


X"(x) = wX(2). (5) 
The boundary conditions (2) imply 
X(0) = X(L) =0. (6) 


For jz > 0 the only function X satisfying (5) and (6) is 
the zero function, leading to the trivial solution U = 0 of 
the original partial differential equation. For 4 < 0a 
non-trivial solution is possible: 

X(x) = Fsinhz, 
where p = —h? (h > 0) and y is chosen so that 
sinhL = 0, ie. 


h 


TL 
where n is a positive integer. With this choice of x, when 
we have solved the equation 

Tr =Aur 
implied by (4), the resulting solution U(x, t) = X(x)T(t) 
takes the form 


U, (x,t) = sin = 


where n = 1,2,3,.... 


= ( nret 


. umet 
An cos FE + sin “Et 


2. A more general solution of the partial differential 
equation is (by Theorem 1): 


: 
nmr niet . nmet 
Dea (40 cos "+ By sin“). (7) 


If the solution is required to satisfy initial conditions of 
the form 
U(x,0) = (x) 
= initial position of point xr on spring, (8a) 
w(z) 
nitial velocity of point xx on spring, (8b) 


U(z,t) = 


ou 
9) 


then the coefficients A, and B, can be found from the 
Fourier sine series expansions for the initial position and 
velocity functions ¢(x) and w(x) defined in Equations 
(8a) and (8b). 


Section 4: Transverse vibrations of a stretched 
string 


1. Transverse vibrations of a heavy spring or string can 
also be modelled by the wave equation, though the value 
of the wave speed c is not the same as for longitudinal 
vibrations. 


2. A boundary condition with the property that if two 
functions U; and U3 satisfy it then so do all linear 
combinations c)U; + ¢2U% is said to be homogeneous; 
otherwise it is non-homogeneous. 
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3. The wave equation (1) with the non-homogeneous set 
of boundary conditions 

U(0,t) = 0, (9a) 

U(L,t) = Rsinut, (9b) 
with w a constant, can be used to model transverse 
vibrations of a stretched string, one end of which is fixed 
and the other constrained to move sinusoidally, 
For this problem, separation of variables gives the 
particular solution 
Rsin(we/c) 
FinwL/e) St: 
The general solution is 

U(x, t) = Uo(x,t) + V(z,t), 
where V(x, t) stands for the right-hand side of Equation 
(7). If initial conditions are given as well, the constants 
Ap, By can be determined using Fourier series, as in 
Section 3. 


Uo(x,t) = 
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9 Index of definitions 


A reference in this index beginning with the letter ‘p’ is to a particular page of this 


Handbook. 


A reference consisting of only numbers refers to the unit-by-unit summaries in 
Section 8 of this Handbook. Normally, the reference consists of two or three numbers, 
of which the first describes the unit, the second the section, and the third (where 
applicable) the numbered item, for those sections which are split into items. 

(For example, 29.3.4 means Item 4 of Section 3 of the entry for Unit 29.) 


A number on its own refers to the introductory part of the entry. (For example, 6 
means the Introduction of the entry for Unit 6.) 


absolute error 

absolute stability of numerical method 

absolute zero (of temperature) 

absolutely ill-conditioned 

absolutely well-conditioned 

AC — B? criterion 

acceleration (one-dimensional) 

acceleration (three-dimensional) 

acceleration vector 

addition of matrices 

addition of vectors 

air resistance 

amplitude 

angular frequency 

angular momentum 

angular velocity 

Argand diagram 

argument of complex number 

associated homogeneous equation of 
differential equation 

associativity 

augmented matrix 

auxiliary equation of differential equation 

auxiliary equation of recurrence relation 

axioms of algebra 


back substitution 
backward recurrence 
binomial theorem 
boundary condition 


Cartesian component of vector 
Cartesian form of complex number 
Cartesian unit vectors 

catenary 

central force 

centre of mass 

centre of mass decomposition theorems 
chain rule for ordinary differentiation 
chain rule for partial derivatives 
characteristic equation of matrix 
circular motion 

closure 

codomain 

coefficient of kinetic friction 
coefficient of restitution 

coefficient of static friction 
coefficients of differential equation 
coefficients of recurrence relation 
collision 

column matrix 


29.1.2, 29.1.4 
28.3.2, 28.4.6 
5.2 
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21.1.2 


column vector 

combined heat transfer coefficient 
commutativity 

complementary function 

completing the square (for a quadratic) 
complex conjugate 

complex number 

component of force 

component of vector 

composite form for integration method 
composite function 

composition of functions 

conditioning of Gaussian elimination 
conduction 

conservation of angular momentum 
conservation of linear momentum 
conservation of (mechanical) energy 


conservative field 

conservative force 

consistent numerical method 
constant acceleration (motion with) 
constant catch strategy 

constant effort strategy 

constant function 

constant velocity (motion with) 
constant-coefficient differential equation 
constant-coefficient recurrence relation 
constraint 

continuous function 

contour 

contour curve 

contour surface 

convection 

convective heat transfer coefficient 
convergent numerical method 
critical damping 

cross product 

cubic (polynomial) 

curl 

cylindrical polar coordinates 


damped harmonic oscillator 
damped oscillations 
damped vibrations 
damping 

damping constant 

damping ratio 

dashpot 

dashpot constant 

dashpot force 


datum of potential energy 

de Moivre's Theorem 

decomposition of matrix 

decomposition theorems 

degeneracy of linear programming 
problem 

degenerate vibrating system 

degree (of a polynomial) 

degrees of freedom 

dependent variable 

derivative 

derivative of vector function 

derived function 

determinant 

diagonal matrix 

diagonal of matrix 

difference of two squares 

differential vector operator 

differentiating partially 

differentiation 

direct iteration 

direction field of differential equation 

direction of vector 

displacement ratio 

displacement vector 

distinct eigenvalues 

distributivity 

domain 

dot product 


eccentricity (of ellipse) 
effective potential energy 
eigenvalue 

eigenvector 

elastic collision 

element of matrix 

ellipse 

empirical model 

equality of matrices 
error bound 

error function 

essential row interchange 
Buler’s formula for e“@ 
Euler's integration formula 
Euler’s method 


Euler-trapezoidal method 

even function 

explicit numerical method 

exploitation of populations 

exponential form of complex number 

exponential function of complex 
number 

exponential model for population 

external force 

external torque 

extrapolation 

extremum 


factorial 
factorizati 
feasible region 

field line 

first-order differential equation 
first-order partial derivative 
first-order recurrence relation 


20.1.2 

18.1.6, 18.4.3 
18.1.5 

9.3.1 

54 

18.4.4 

2.2, 6.5, 
19.1.2, 22.2.2 


force of friction 
force of gravity 


force (one-dimensional) 

force (three-dimensional) 

force vector 

forced oscillations 

forced vibrations 

forward recurrence 

Fourier series 

Fourier’s law for conduction 
Fourier’s Theorem 

frequency 

frequency response 

friction 

frictional force 

frictional reaction 

function of a function rule 
function of many variables 
function of one variable 
fundamental angular frequency 
fundamental theorem of algebra 
fundamental theorem of calculus 


Gaussian elimination 


general solution of differential equation 
general solution of recurrence relation 


global error 

grad 

gradient of scalar function 
graph of function 
graphical method 
gravitational constant 
gravitational force 


gravitational potential energy 
gravity (force of) 


heat transfer 

heat transfer coefficient 
homogeneous boundary condition 
homogeneous differential equation 


homogeneous linear algebraic equations 


homogeneous recurrence relation 
Hooke’s law 

hyperbolic cosine function 
hyperbolic sine function 


identity element 

identity matrix 

image of a function 

imaginary (complex) number 

imaginary part of complex number 

implicit differentiation 

implicit numerical method 

independent variable 

induced instability of Gaussian 
elimination 

inelastic collision 

inhomogeneous differential equation 

inhomogeneous recurrence relation 


initial condition for differential equation 
initial condition for recurrence relation 


input-output principle 
integrating factor 
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28.4.2, 


4.3.2, 28.4.2, 


28.4.2, 


12.4.1, 12.4.4 
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integration 
integration by parts 
integration by substitution 
internal force 

inter-particle force 
interpolating polynomial 
interpolation 

interval of absolute stability 
inverse (element) 

inverse iteration 

inverse matrix 

isotropic force 


Kepler's first law 

Kepler's laws 

Kepler's second law 

Kepler’s third law 

kinetic energy 

kinetic energy of rotating rigid body 


law of gravitation 

length of curve 

limit of a function 

line integral 

linear combination of functions 
linear combination of matrix rows 
linear damping, 

linear differential equation 
linear function 

linear momentum 

linear polynomial 

linear programming problem 
linear recurrence relation 
linear system of differential equations 
linear transformation 

linearity rule for integration 
linearly dependent matrix rows 
linearly independent functions 
local maximum 

jocal minimum 

local truncation error 

logistic curve 

logistic equation 

logistic model 

longitudinal vibrations 

lower triangular matrix 

LR method for eigenvalues 

LU decomposition 
lumped-parameter model 


magnitude of vector 

main diagonal of matrix 
matrix 

matrix addition 

matrix element 

matrix equality 

matrix multiplication 
matrix of coefficients 
matrix subtraction 
maximum 

Mean Value Theorem 
mechanical energy 
minimum 

Mixed Derivative Theorem 
modelling (mathematical) 


30.2.1 


modified inverse iteration method 
modulus of complex number 
modulus of elasticity 
modulus of vector 

moment of force 

moment of inertia 
momentum (angular) 
momentum (linear) 
mortgage problem 
multiplication of matrices 
multipliers 


natural length of spring 

nested multiplication for polynomials 

Newton-Cotes formula 

Newton-Raphson method 

Newton's first law 

Newton's law of restitution 

Newton's second law (one-dimensional) 

Newton's second law (three-dimensional) 

Newton's third law 

Newton's universal law of gravitation 

non-homogeneous differential equation 

non-homogeneous recurrence relation 

non-singular matrix 

normal form of system of differential 
equations 

normal mode 

normal mode angular frequency 

normal reaction 

normal vector 


objective function 10.1 
odd function 31.3.1 
one-step numerical method 19.3.1 
optimal vertex 10.2 
optimum 10.1 
orbit 30.3.4 
order of differential equation 241.1, 6 
order of recurrence relation 1.1.2 
parallel axes theorem 29.2.5 
parametric equation 26.3.2 


partial derivative 

partial differential equation 

partial differentiation 

partial fractions 

partial pivoting 

particle 

particular solution of differential equation 
particular solution of recurrence relation 
path of integration 

perfect dashpot 

perfect spring 

perfect square (polynomial) 

period 


periodic function 
phase (angle) 

phase shift 

phasor 

piecewise continuous 

pivot for Gaussian elimination 


pivot row 3. 
plane polar unit vectors 28.4.4, 30.1.3 
polar coordinates p. 10 


polar form of complex number 

polar unit vectors 

polynomial 

position (one-dimensional) 

position (three-dimensional) 

position vector 

potential energy 

power 

precession 

predictor-corrector method 

principal term in truncation error 

principal value of angle 

principal value range of inverse 
trigonometric function 

principle of superposition 

product rule for differentiation 

projection of vector 

proportionate birth rate 

proportionate death rate 

proportionate growth rate of 
population 


quadratic (polynomial) 
quadratic Taylor approximation 
quotient rule for differentiation 


radial component 

radiation 

range of projectile 

rate of precession 

rational function 

reaction 

real (complex) number 

real part of complex number 
recurrence relation 

region of integration 
relative error 

relative response 

relatively ill-conditioned 
relatively well-conditioned 
remainder for approximation 
repeated eigenvalue 
resonance 

restitution 

right-hand rule 

right-hand screw rule 
right-handed coordinate system 
rigid body 

rocket equation 

rolling condition 

root (of a polynomial) 

root of complex number 
rotation of rigid body 

row matrix 

row vector 

rule (of a function) 


saddle point 
scalar 

scalar field 

scalar line integral 
scalar product 
scale factor 
scaling of vector 
Schaefer model 


5.3 
28.4.4, 30.1.3 


LL 
27.1.1, 27.2.1 


second derivative 

second-order differential equation 

second-order partial derivative 

second-order recurrence relation 

section function 

semi-latus rectum 

semi-major axis 

semi-minor axis 

separation of variables for first-order 
differential equation 

separation of variables for partial 
differential equation 

SI units 

simple harmonic motion 

simplex method 

Simpson's integration formula 

Simpson's method 


simultaneous linear algebraic equations 


singular matrix 

sinusoidal function 

size of matrix 

specific energy capacity 
specific heat 

speed 

spherical polar coordinates 
spring force 

spring potential energy 
spurious solution 

square matrix 

stability of Gaussian elimination 
stability of numerical method 
static friction 

static rigid body 

stationary point 

steady-state vibration 

step length 

step size 

stiffness of spring 

strong damping 

substitution rule for integration 
subtraction of matrices 
subtraction of vectors 

sum rule for differentiation 
summation notation 
superposition 

surface 

surface integral 

symmetric matrix 


tangent approximation 
tangent plane 
Taylor approximation 


Taylor polynomial 


Taylor series 

Taylor series method 

Taylor's Theorem 

temperature gradient 

tension of spring 

terminal speed 

thermal conductivity 

third derivative 

thrust of spring 

torque law about centre of mass 


p.13 
6 

25.3.1 

1.2, 5.5 
25.1.3 

30.2.2, 30.3.4 
30.2.2 

30.2.2 


15.3.4 


25.4.3 


p. 15 
25.2.4 
p. 15, 18,1.3, 
25.2.2, 25.3.4 
p. 15, 18.1.3, 
25.2.1, 25.3.3 
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torque law about fixed point 


torque (of force) 

trajectory of projectile 

transfer coefficient 

transient 

transpose of matrix 

transverse component 

transverse vibrations 

trapezoidal integration formula 

trapezoidal method 

triangle inequality 

triangular matrix 

trigonometric functions 

trigonometric identities 

triple scalar product 

trivial constraint 

trivial solution of homogeneous 
linear equations 

truncation error 

turning points of motion 

two-step numerical method 


U value 

undamped angular frequency 
uniform circular motion 
uniform motion 

unit matrix 

unit vector 

universal law of gravitation 
upper triangular matrix 


vector 

vector addition 

vector differential operator 
vector field 

vector function 

vector product 

vector subtraction 

velocity (one-dimensional) 
velocity (three-dimensional) 
velocity vector 

vibration absorber 

volume integral 


wave equation 
wave speed 
weak damping 


yield of exploitation of populations 


zero function 
zero matrix 
zero vector 


29.1.3, 
29.1.6, 29.2.3 
28.4.1, 29.1.5 
15.4.1 

12.4.1, 12.4.4 
6.4 

20.2.4 

28.3.1 

32.2.3, 32.4.1 
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10 Notation 


The following symbols have standard meanings throughout the course. 


i square root of —1 (see 5.1) 

g magnitude of acceleration due to gravity 
(approximately 9.81 ms~) (see 4.3.2) 

sin, cos, tan p.8 

sec, cosec, cot p.8 

aresin, arccos, etc. p, 9 

grad, curl see Unit 26 

lim p. ll 

log, x the logarithm of x to the base e = 2.71828... 

> approximately equals 

! factorial 

a/be a+(bxc) 

ve the positive square root of x 

alr nth root of number x 

as 1/2” 

|a| modulus of real number x (= Vz?) 

lel modulus of complex number = 

z complex conjugate of complex number = 

|v| or v modulus or length of vector v 

v unit vector along v 

w(x), f(x) see p. 11 

tid derived function of function f (p. 13) 

y” derived function of f’ (p. 13) 

wis) nth derivative of f (p. 13) 

ay first derivative, i.e. f’(x) where y = f(x) 

& &. ie. f’(t) where x = f(t) 

2 oe a 

az ef ie. f”(t) where x = f(t) 

v ie. s"(2) where y = f(x) 

Byn a wy Bee, el ay 

mo y" second derivative, i.e. i: (#) 

ou or y™) nth derivative 

/ f(x) dx indefinite integral of function f(x) 

b 
f(x) dx definite integral of function f(x) from z =a tor =6 
a 

€ belongs to 

{a, 6] the set of all x such thata <2 <b 

(a,b) the set of all x such that a<2<b 

= summation (see p. 5) 

Or a-axis of coordinate system 


(x, y)-plane plane containing «- and y-axes 
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